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Abstract 


>*The  binary  error  probability  for  integrate-and-dump 
receivers  is  developed.  The  combination  of  the  optical 
modulator#  optical  channel#  and  optical  detector  is  charac¬ 
terized  as  a  discrete  memoryless  channel  (DMC) .  Maximum  a 
posteriori  (MAP)  decision  criteria  are  developed  to  enable 
optimal  post-detection  processing.  The  effect  of  quantizing 
the  output  is  addressed  for  both  ideal  unity  gain  detector 
and  avalanche  detector  optical  systems.  Multilevel  Amplitude 
Keying  is  used  to  calculate  specific  probability  of  error 
values#  and  the  channel  cutoff  rate  is  used  to  investigate 
probability  of  error  bounds. 


x 


INVESTIGATION  OF  BINARY  ERROR  PROBABILITY 
FOR  AN  INTEGRATE-AND-DUMP  RECEIVER 
AS  OUTPUT  QUANTIZATION  INCREASES 


I 

i 

i 

1  *1  Introduction 


y 

This  thesis  effort  focuses  on  direct  detection  optical 
communications  systems  which  use  receivers  classified  as 
photon  counters,  or  integrate-and-dump  receivers.  These 
systems  determine  the  received  message  estimate  based  upon 
the  number  of  photons,  or  counts,  accumulated  in  the  speci¬ 
fied  reception  time  period,  T.  Figure  1  shows  a  typical 
block  diagram  of  an  optical  communication  system. 

This  report  centers  on  the  examination ■ of  the  optical 
channel,  the  optical  detector,  and  the  post-detection  pro¬ 
cessor,  or  decision  maker.  The  optical  channel  and  detector 
are  characterized,  or  modeled,  and  the  post-detection  pro¬ 
cessing  is  then  optimized  for  the  channel  model  to  provide 
minimum  probability  of  error  for  the  system  as  a  unit.  The 
channel  and  detector  are  analyzed  together  and  are  modeled 
as  a  discrete  memoryless  channel  ( DMC )  with  characteristic 
crossover  or  transition  probabilities.  Using  the  channel 
model,  the  optimum  post-detection  processor  is  developed 
using  maximum  a  posteriori  (HAP)  decision  criteria  and  the 
resulting  system  is  evaluated  by  determining  the  probability 
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Figure  1.  Optical  Communication  System  Block  Diagram 


of  error.  The  effect  of  quantizing  the  output  of  the  detec¬ 
tor,  or  dividing  the  time  period  T  into  smaller  segments,  is 
examined. 

Throughout  this  work,  a  wideband  direct  detection 
receiver  has  been  assumed  to  allow  the  signals  to  be  viewed 
as  point  processes.  This  assumption  essentially  limits  all 
disturbances  to  those  associated  with  the  conversion  of  the 
optical  field  to  photons  (Ref  2:113).  A  binary  signal  input 
field  is  assumed,  however,  the  procedures  used  are  easily 
generalized  to  an  M-ary  input  field. 

The  concepts  used  to  analyze  the  basic  optical  communi¬ 
cation  system  of  interest  are  developed  in  Chapter  II.  The 
general  procedure  used  to  characterize  the  channel/detector 
combination  is  presented,  and  the  selection  method  for  the 
optimum  post-detection  processor  for  unquantized  output  is 
examined.  Chapter  III  introduces  the  idea  of  quantized 


output  and  applies  the  concepts  developed  in  Chapter  II  to 
communications  systems  employing  increasingly  finer  output 
quantization.  The  cutoff  rate  approach  to  communications 
systems  is  presented  in  Chapter  IV.  The  meaning  of  the 
cutoff  rate  is  examined,  and  probability  of  error  bounds  are 
presented.  Chapter  V’ introduces  a  random  gain  factor  to  the 
output  of  the  optical  detector,  modeling  the  effect  of 
including  an  avalanche  photodiode  in  the  optical  detection 
system.  The  effect  of  quantizing  the  output  of  the  avalanche 
detector  is  examined  by  calculating  probability  of  error 
information.  The  final  chapter  summarizes  the  major  results 
of  this  research  and  discusses  areas  open  for  future  research 


II  Background 


tical  Communication  System  Model 


A  typical  optical  communication  system  is  shown  m 
Figure  1.  Some  type  of  information  is  to  be  transmitted 
from  a  source  to  a  distant  destination  by  propagation  of  a 
modulated  light  signal.  The  information  signal  will  vary  a 
specific  property  of  the  transmitted  light  such  as  amplitude 
frequency,  or  phase.  The  optical  field  detected  at  the 
receiver  depends  directly  on  the  field  that  was  transmitted, 
the  effects  introduced  by  the  channel,  and  background 
radiation.  The  function  of  the  receiver  is  to  process  the 
received  optical  field  so  as  to  reproduce  an  accurate  esti¬ 
mate  of  the  information  signal. 


Integrate-and-Dump  Receivers 


The  development  of  the  theory  applicable  to  direct 
detection  receivers  has  been  well  documented  (Refs  2,  6,  8) 
Given  a  field  incident  on  the  detector  surface,  Ud  U,t>, 
the  rate  function,  X(t),  is  defined  as: 

A(t)  =  Cyhf0)  /  fLfd(r,t)i2  dr  (1) 

Ao 

where 


detectors  quantum  efficiency 
energy  of  a  photon 


detectors  active  surface  area 


Figure  2.  Integrate-and-Dump  Receiver  Block  Diagram 


This  rate  function,  X(t),  gives  the  number  of  photons  inci¬ 
dent  on  the  detector,  per  time  period,  for  a  given  optical 
field.  In  optical  communications  systems  which  use  photon 
sensitive  detectors,  the  emission  of  the  photoelectrons 
corresponds  to  a  point  process  obeying  Poisson  statistics 
(Refs  1,  18). 

The  block  diagram  of  an  integrate-and-dump  receiver 
shown  in  Figure  2  gives  an  overall  view  of  how  the  counting 
of  photons  is  implemented.  By  integrating  the  output  of  the 
detector,  the  rate  function  X(t),  over  the  time  period  T,  a 
total  number  of  counts  for  T,  N(T) ,  is  obtained.  This 
number,  N(T),  reflects  the  number  of  photons  emitted  by  the 
detector  in  time  T.  Thus,  as  different  optical  fields, 
l — r(r,t),  will  be  associated  with  each  of  the  i  possible 
source  messages,  there  will  also  be  different  rates,  X^(t), 
and  counts,  N^(T),  associated  with  each  source  message. 
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Channel  Characterization 

As  demonstrated  by  Snyder  and  Rhodes  (Ref  17:327),  the 
combination  of  the  optical  modulator,  optical  channel,  and 
the  optical  detector,  when  considered  as  a  unit,  forms  a 
discrete  channel.  For  the  purpose  of  this  study,  the  chan¬ 
nel  will  be  assumed  to  be  memoryless,  that  is,  the  channel 
output  at  any  given  time  is  a  function  of  the  channel  input 
at  that  time  and  is  independent  of  all  previous  inputs.  A 
discrete  memoryless  channel  ( DMC )  can  be  completely  specified 
by  the  set  of  conditional  probabilities  which  relate  the 
probability  of  each  output  to  the  various  inputs  (Ref  7:6). 

The  assumption  of  a  binary  input  signal  set  leads  to 
the  use  of  a  binary  input  channel  shown  in  Figure  3.  By 
combining  the  modulator,  channel,  and  detector  to  form  a  DMC 
and  determining  the  probabilities  pQ  and  p^  of  Figure  3,  it 
becomes  possible  to  determine  the  optimum  post-detection 
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processing  that  will  be  required  to  produce  minimum  proba¬ 
bility  of  error  estimates.  / 

t 

C  *- imu  Post-Detection  Processing  j 

T  e  optimum  post-detection  processing  is  that  process¬ 
ing  of  the  channel  output  required  to  provide  an  estimate  of 
the  input  information  signal  while  maintaining  the  system 
probability  of  error  at  a  minimum.  Essentially,  the  post¬ 
detection  processor  is  a  decision  maker  that  determines  what 
the  input  information  signal  was,  based  on  the  output  of  the 
channel . 

The  optimum  post-detection  processing  reduces  to  simple 
photon  counting  for  the  following  cases: 

1  -  When  Xi(t)  is  a  constant  for  t  within  the  specified 

time  period  of  interest,  t  e[o,T]; 

2  -  When  X.(t)  is  the  same  for  all  signals,  but  the 

signals  are  spatially  separated;  or, 

3  -  When  X^(t)  are  nonoverlapping,  or  orthogonal,  in 

time. 

The  first  case  is  commonly  called  Multilevel  Amplitude 
Keying,  while  the  third  case  is  Pulse  Position  Modulation 
(Ref  2:171).  This  report  will  give  examples  for  Multilevel 
Amplitude  Keying.  By  assuming  Multilevel  Amplitude  Keying, 
the  dark  current,  which  is  a  homogeneous  Poisson  process 
with  rate  X^,  can  be  included  within  the  rates  X^^  because 
the  dark  current  will  have  the  same  effect  on  any  of  the 
received  signals.  That  is,  the  rate  XD,  when  integrated 
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over  time  period  T,  will  add  a  certain  number  of  counts 
independent  of  which  source  message  was  transmitted. 

To  ensure  optimum  post-detection  processing,  that  is, 
minimum  probability  of  error  for  the  system,  the  maximum  a 
posteriori  (MAP)  decision  criterion  will  be  used.  The  MAP 
decision  criterion  states  that  given  an  observation,  Y, 
select  decision  one  if  M^  is  more  probable  than  Mg  to  have 
been  the  input  information  signal  which  produced  Y  (Ref 
13:42).  Mathematically,  the  decision  criterion  can  be 
written: 


Pr{M. ! Y }  °1 


By  using  Bayes  rule,  Eq  (2)  can  be  rewritten  as: 


PrlYjMj}  ®1  Pr{M0) 


This  introduces  the  concept  of  decision  regions.  Whenever 
the  left-hand  side  of  Eq  (3),  commonly  called  the  likelihood 
ratio,  is  greater  than  the  threshold,  given  by  the  right- 
hand  side  of  Eq  (3),  the  decision  is  made  that  M^  was  the 
information  transmitted.  Thus,  the  threshold  divides  space 
into  two  decision  regions:  one  where  the  likelihood  ratio 
is  greater  than  the  threshold  and  another  where  it  is  less. 

Returning  to  the  ideal  photon  detector,  the  optimum 
post-detection  processor  will  decide  that  information 
message  one  was  transmitted  whenever  the  probability  that 


Ideal 


t=T 


Figure  4.  Binary  Input,  Unquantized  Output  System 


the  number  of  counts  obtained,  N ( T ) ,  was  more  likely  given 
message  one  than  message  zero. 

Determining  Probability  of  Error 

Once  the  decision  regions  have  been  determined,  the 
system  probability  of  error  can  be  calculated.  The  total 
probability  of  error  is  simply  the  probability  of  deciding 
message  zero  when  message  one  was  sent  plus  the  probability 
of  deciding  message  one  when  message  zero  was  sent.  Mathe¬ 
matically: 

P(E)  =  Pr{D0,M1>  +  Pr{DirM0)  (4) 

and,  using  conditional  probabilities: 

P(E)  =  Pr {Dq | M1)Pr {M1J  +  Pr{D1|M0)Pr{M0)  (5) 

Binary  Input ,  Unquantized  Binary  Output 

The  system  shown  in  Figure  4  is  a  binary  input,  unquan 
tized  output  optical  communication  system  receiver.  A 
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binary  input  alphabet  of  Mg  and  produce  optical  fields 
H  and  h  incident  to  the  ideal  detector.  The  output  of 
the  ideal  detector  is  the  rate  function,  either  X q ( t )  or 
A^(t).  The  rate  function  is  integrated  over  the  time  period 
T,  producing  a  count,  I J ( T ) .  This  count,  or  the  number  of 
photoelectrons  produced,  is  input  to  the  decision  maker, 
which  using  the  MAP  decision  criterion,  will  determine  which 
message.  Mg  or  M^,  was  transmitted. 

The  characterization  of  the  channel  is  the  first  step 
to  be  examined.  As  previously  discussed,  for  an  ideal 
detector,  the  count  N(T)  is  distributed  as  a  Poisson 
process.  Thus: 


p"  e‘* 


Pr{N(T)  -  n}  =  n! 


where 


V  -Q/  X(a)  da 


Equation  (7)  determines  the  parameter  y,  the  average  number 
of  counts  in  time  T. 

For  the  unquantized  detector  output,  the  parameter  V  is 
determined  by  integrating  X(t)  over  the  time  period  T. 

Since  each  input  signal,  Mg  and  M^,  will  have  a  correspond¬ 
ing  rate,  Ag(t)  and  A^t),  two  y  parameters,  yQ  and  y^,  will 
be  defined.  Applying  these  values  for  y  into  Eq  (6)  gives: 

V0n  e-yo 


Pr{N(T)  *  n|M0) 


(8) 


Pr{N(T)  =  n | Mx }  = 


n  -y, 
e 


Now  that  the  probability  of  receiving  N ( T )  =  n  counts  is 

known  for  any  n,  given  either  Mg  or  M^  was  transmitted, 

« 

Eq  (3)  can  be  used  to  define  the  decision  regions.  Applying 
the  MAP  decision  rule  of  Eq  (3)  to  the  case  of  unquantized 
outputs  given  by  Eqs  (8)  and  (9)  results  in: 


Since  nf  the  number  of  counts,  is  the  variable  of  interest, 
the  decision  regions  should  be  defined  as  functions  of  n. 

To  determine  the  value  for  n  of  interest,  Eq  (10)  is  manipu 
lated  as  follows: 


(^)( 


ni 

y0n  e’*0 


Pr{M0) 
Pr{M, } 


u”  e'yl  ?1  Pr(M„} 


t 


i 


,v 

»*' 
,  w* 


4 


1 


,  /  Pr{Mo}  \ 

eCvtl  ‘  y0)  (  Pr{M1J  J 


(14) 


(15) 


(16) 


for  all  >  Mg  .  Thus,  Eq  (16)  gives  the  decision 
regions  based  on  n,  the  number  of  counts.  The  threshold 
value  is  given  by  the  right-hand  side  of  Eq  (16): 


Hn)  + 


In 


(  pr<M0>\ 

\phm p~y 


(17) 


If  N(T)  is  greater  than  Nq,  the  decision  is  made  that  M1  was 
transmitted.  Alternatively,  if  N(T)  is  less  than  NQ,  the 
decision  is  made  that  Mg  was  transmitted. 

The  count,  N(T),  will  always  be  an  integer  value,  while 
the  threshold,  Ng,  is  continuous  and  may  take  on  any  value. 
For  cases  where  Ng  is  not  an  integer,  the  decision  regions 
are  defined  by  deciding  Mg  for  all  N(T)  which  take  on 
integer  values  less  than  Ng  and  deciding  M^  for  all  N(T) 
which  take  on  integer  values  greater  than  Ng.  If  NQ  is  an 
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integer  value,  the  case  of  N(T)  =  Ng  must  also  be  consid¬ 
ered  since  the  probability  of  N (T)  *  Ng  is  no  longer  zero, 
as  is  the  case  when  Ng  is  not  an  integer.  The  question  of 
which  decision  is  to  be  made  when  on  the  boundary  of  the 
decision  regions  is  answered  as  follows — Choose  for  8 
percent  of  the  time  when  N(T)  =  Ng  ,  where  8  is  given  by: 


a  -  Pr{N(T)  >  Ng|Mg} 
Pr {N (T)  =  N0|M0) 


(18) 


and  a  is  a  parameter  chosen  by  the  user,  satisfying: 

Pr  (N (T)  >  N0|MQ}  <  a  <  Pr{N(T)  >  Ng|Mg}  (19) 

(Ref  2:125D-125I) .  While  may  be  chosen  arbitrarily  when 
N(T)  *  Ng  ,  the  procedure  outlined  above  allows  the  user  to 
determine  what  performance,  or  probability  Of  false  alarm, 
is  acceptable  by  varying  the  parameter  a. 

Thus,  the  optimum  post-detection  processor  makes  the 
decision  between  Mg  and  M^  by  comparing  the  received  count, 
N(T),  to  the  threshold  value,  Ng,  and  choice  of  MQ  or 
depends  on  whether  N  (T )  is  less  than,  greater  than,  or  equal 

tO  Ng. 

The  probability  of  error  for  the  binary  system  is  given 
by  Eq  (5).  The  prior  probabilities  of  Mg  and  are  assumed 
to  be  known,  and  the  probabilities  of  making  an  incorrect 
decision  are  given  by: 


Pr{Dj | Mg }  =  Pr{N(T)  >  Ng|Mg} 


(20) 


TT1 


Pr{Dj |Mq)  * 


0 


n! 


(21) 


and 

Pr{DQ  |*M1 }  =  Pr {N (T)  <  NgjMj}  (22) 

or 

LNoJ  y^  e"yl 

Pr{D0lMl}  =  nlQ  -  (23) 

Therefore,  all  of  the  terms  of  Eq  (5)  are  known  and  the 
probability  of  error  can  be  calculated. 

Summary 

The  calculation  of  binary  error  probability  for  an 
integrate-and-dump  receiver  is  accomplished  in  a  straight¬ 
forward  manner.  The  first  step  is  to  characterize  or  model 
the  channel.  Then,  using  the  channel  model,  the  optimum 
post-detection  processor  for  that  channel  must  be  determined 
Once  these  steps  have  been  completed,  the  probability  of 
error  is  a  direct  result. 

For  the  unquantized  output  case,  the  MAP  decision  rule 
is  determined  directly  from  the  known  values  for  the  prior 
probabilities  of  Mg  and  M^  and  the  average  counts  received 
for  each  message,  Ug  and  y^.  Table  I  lists  calculated 
values  of  the  probability  of  error  for  various  values  of 


prior  probabilities  and  average  counts  per  message.  Appendix 
A,  Figure  15  is  a  listing  of  the  BASIC  language  computer 
program  used  to  calculate  the  data.  Appendix  B,  Table  II  is 
a  listing  of  calculated  data  for  both  quantized  and  unquantized 
output,  ideal  detector  systems. 
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TABLE  I 


Probability  of 

Error  - 

Unquantized 

Output 

Pr^} 

VO. 

ILL 

Ifl. 

P  ( E ) 

0.9 

1 

4 

0.58 

0.0797 

0.9  * 

1 

5 

1.12 

0.0628 

0.7 

1 

4 

1.55 

0.1434 

0.7 

1 

5 

1.96 

0.1076 

0.5 

1 

4 

2.16 

0.1592 

0.5 

1 

5 

2.49 

0.1025 

0.3 

1 

4 

2.78 

0.1276 

0.3 

1 

5 

3.01 

0.0928 

0.1 

1 

4 

3.74 

0.0604 

0.1 


1 


5 


3.85 


0.0436 


Ill  Increasing  Output  Quantization 


Introduction 

Having  already  developed  the  procedure  for  determining 
the  optimum  post-detection  processing  scheme  and  the  resul¬ 
tant  probability  of  error  for  unquantized  output,  ideal 
detector  receivers,  this  section  introduces  output  quantiza¬ 
tion  into  the  model.  The  characterization  of  the  channel 
will  be  re-examined,  the  optimum  post-detection  processor 
will  be  redefined,  and  the  resulting  probability  of  error 
will  be  calculated  for  varying  levels  of  output  quantization 

Channel  Characterization 

For  the  unquantized  output  communication  system,  the 
channel  model  resulted  in  a  count  of  photoelectrons,  N(T). 

As  output  quantization  is  introduced,  the  time  period  T  is 
divided  into  shorter  length  segments.  By  doing  this,  the 
count  N(T)  is  also  broken  up  into  several  counts,  one  count 
for  each  segment  in  T.  For  example,  if  T  is  segmented  into 
two  segments,  each  of  (T/2)  duration,  then  there  will  now  be 
two  counts  associated  with  T,  one  for  each  (T/2)  time 
interval.  In  general,  if  T  is  segmented  into  J  intervals, 
each  of  duration  (T/J),  then  there  will  be  J  individual 
counts  to  be  considered  in  the  decision  process.  These  J 
individual  counts  will  be  symbolized  by  N^,  for  i  *  1  to 
J,  and  J  is  the  quantization  level  of  the  output.  It  is 


important  to  note  that  the  unquantized  system  is  obtained  by 
setting  J  equal  to  one. 

By  quantizing  the  output,  the  decision  will  now  be 
based  on  J  bits  of  information,  the  counts,  rather  than 
the  single  count,  N(T).  By  using  a  MAP  decision  rule  on 
each  count,  ,  and  deciding  Mq  or  for  each  count,  a  J 
bit  long  string  of  zeros  and  ones  can  be  developed.  It  is 
this  J  bit  long  string  that  will  be  the  output  of  the 
channel  and  the  input  of  the  decision  maker  for  deciding 
whether  Mq  or  M1  was  transmitted,  and  the  optimum  post¬ 
detection  processor  must  make  the  decision  with  minimum 
probability  of  error. 

The  channel  characterization  for  the  quantized  output 
system  requires  the  determination  of  the  transition  probabil¬ 
ities  associated  with  each  of  the  2J  possible  channel  out¬ 
puts  of  J  bit  strings  of  zeros  and  ones.  Figure  5  shows  a 
possible  channel  characterization  for  J  *  2. 

As  was  the  case  for  unquantized  output,  the  counts 
for  i  *  1  to  J,  obey  Poisson  statistics  as  given  by  Eq  (6). 
Due  to  the  quantization,  however,  Eq  (7)  is  modified  to: 

(n+l)T/J 

V-  T  _  *  /  X.(a)da  for  n=0  to  (J-l) 

1,J’n  nT/J  1  (24) 

because  X^(t)  is  no  longer  integrated  over  the  entire  period 
T,  but  only  for  intervals  of  (T/J)  duration.  By  assuming 
Multilevel  Amplitude  Keying,  Eq  (24)  reduces  to: 


Figure  5.  Channel  Characterization  for  J  -  2 
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Now,  the  MAP  decision  rule  of  Eq  (3)  can  be  applied  to  each 


of  the  J  segments  giving: 


Pr{Ni  =  nlMj} 
Pr{Ni  =  n | MQ } 


>1  Pr{M0) 
<  TrlMTT 

n  J- 


(26 


Pr{N.  -  n|M, }  « 


CVJ)n 


where 


(28) 


n  -(P0/J) 

Pr{N-  =  n|Mn)  =  ^0/J)  e _ 

1  u  n! 

Substituting  Eqs  (27)  and  (28)  into  Eq  (26)  and  simplifying, 
as  was  done  previously  in  Eqs  (11)  through  (16),  the  MAP 
decision  rule  for  each  of  the  J  segments  becomes: 


Thus,  for  each  of  the  J  segments  of  T,  if  the  count  for 
a  segment,  N.,  is  greater  than  the  threshold,  N_  T,  set  the 

l  u ,  J 

bit  for  that  segment  to  zero.  For  Ni  =  No  j  '  the  Proce“ 
dure  discussed  in  Chapter  II  using  Eq  (18)  is  used  to  make 
the  decision.  Each  count  for  the  J  segments  is  evaluated 
this  way  until  a  J  bit  long  string  of  zeros  and  ones  has 
been  formed.  This  string  then  represents  the  output  of  the 
quantized  channel. 

To  fully  characterize  the  channel,  the  transition 
probabilities  for  each  of  the  2J  possible  channel  outputs 
are  required.  To  determine  the  transition  probabilities  for 
the  J  bit  long  strings,  the  crossover  probability  for  each 
of  the  J  segments  must  be  known.  That  is,  the  probability 
of  setting  a  segment's  bit  to  a  one  given  Mg  was  transmitted 
and  the  probability  of  setting  a  segment's  bit  to  a  zero 
given  was  transmitted  must  be  determined.  The  procedure 
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to  determine  these  bit  transition  probabilities  is  the  same 
as  was  used  to  calculate  the  system  probability  of  error  in 
Chapter  II.  The  results,  similar  to  Eqs  (21)  and  (23),  are: 


PrtDjlMjj}  = 


77  e 


~T 


n=  FU 


n! 


-  d  -  p0> 


(30) 


and 


PrtDdlMj) 


KJ 


z 

n=0 


n! 


=  (1  -  Px) 


(31) 


These  determine  the  single  bit  transition  probabilities. 

For  an  arbitrary  J,  the  transition  probabilities  associated 
with  a  particular  J  bit  long  string  are  simple  to  calculate. 
The  steps  are  as  follows: 

1  -  Let  the  number  of  ones  in  the  J  bit  string  equal  v, 

then  the  number  of  zeros  in  the  string  is  J  -  v; 

2  -  To  determine  the  transition  probability  of  the 

particular  string,  given  was  transmitted: 


Pr { (J  bit  string)  *  pxV(l  -  P1)J  '  V  (32) 

where  (1  -  p^)  is  given  by  Eq  (31); 

3  -  To  determine  the  transition  probability  of  the 
particular  string,  given  Mn  was  transmitted: 


(33) 


o 


Pr{ (J  bit  string) | MQ >  =  pQJ  V(1  -  pQ)V 

where  (1  -  p^)  is  given  by  Eq  (30). 

At  this  point,  a  model  which  fully  characterizes  the 
channel  has  been  developed.  The  "hard  decision"  channel  of 
the  unquantized  output  system,  where  the  number  of  inputs  to 
the  channel  equals  the  number  of  channel  outputs,  has  become 
a  "soft  decision"  channel  for  the  quantized  output  system, 
where  there  are  a  larger  number  of  possible  channel  outputs 
than  inputs.  Massey  (Ref  11)  and  Lee  (Ref  10)  show  that 
probability  of  error  criterion  is  applicable  only  for  "hard 
decisions*  and  that  a  cutoff  rate  based  criterion  offers 
more  guidance  for  "soft  decision"  system  design.  This  topic 
is  covered  in  Chapter  IV.  The  next  step  to  be  covered  for 
the  quantized  output  system  is  how  to  determine  what  post¬ 
detection  processing  is  required  to  produce  optimum  results. 

Post-Detection  Processing 

Once  the  channel  is  characterized,  the  post-detection 
processor  which  will  take  the  J  bit  string  of  zeros  and  ones 
output  by  the  channel  and  determine  whether  Mg  or  was 
transmitted  with  a  minimum  probability  of  error  must  be 
developed.  Figure  6  is  a  block  diagram  of  this  type  of 
decision  maker. 

To  ensure  minimum  probability  of  error,  the  MAP 
decision  criterion  of  Eq  (3)  will  be  used  once  again.  For 
any  J  bit  string  of  zeros  and  ones  output  by  the  channel. 
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J  bit  lone  strines 


Decision  Maker 


one  string 


2 J  possible  strings 


Figure  6.  Quantized  Output  Post-Detection  Processor 


implementing  the  MAP  decision  rule  of  Eq  (3)  yields: 


Pr{X1X2. . . X j I M1 >  "l  Pr{MQ} 

Pr(X1X2. . .Xj |M0)  <  Pr{M1) 

D  . 


(34) 


where  X^X2  .  .  .  Xj  is  the  J  bit  long  string  of  zeros  and 
ones.  Substituting  Eqs  (32)  and  (33)  into  Eq  (34)  results 


p/d  -  Pj)J  '  V 

PoJ  ‘  VC1  •  Po)V 


Pr{MQ} 

TrTMTT 


(35) 


where  v  is  the  number  of  ones  in  the  J  bit  string.  Similar 
to  the  steps  taken  to  determine  NQ  j  in  Eq  (29),  Eq  (35)  can 
be  simplified,  yielding: 


V 


< 


L 


0 


(36) 


Thus,  the  decision  maker  of  Figure  6  determines  whether 

« 

Mq  or  was  transmitted  by  deciding  if  v,  the  number  of 
ones  in  the  J  bit  string  is  greater  than  or  less  than  Lp. 

For  all  J  bit  strings  where  the  number  of  ones  is  greater 
than  the  threshold,  Lg,  is  decided  to  have  been  transmit¬ 
ted.  If  the  number  of  ones  is  less  than  Lg,  then  MQ  is 
decided.  For  v  =  Lg  ,  the  procedure  discussed  previously 
in  Chapter  II  using  Eq  (18)  is  used  to  make  the  decision. 


Probability  of  Error 

The  system  probability  of  error,  P(E),  is  again  given 
by: 


P(E)  =  Pr {Dj |Mg}Pr {Mg}  +  Pr{DQ |M1}Pr{M1>  (5) 

For  systems  using  quantized  output,  Eq  (5)  reduces  to: 


P(E)  = 


N/j\  v 

Mv]Pl  ci-^ 


J  -  v 


PrtMj} 


J  J 


J- 

v-  M 


_  J  -  v,. 
P0  d 


Po^ 


Pr{M0) 


(37) 


where 
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v/  =  vl(J-v)!  =  binomial  coefficients 


(38) 


Summary 

As  the  output  is  quantized,  the  calculation  of  the 
probability  of  error  becomes  more  complex.  Following  the 
methods  used  for  unquantized  systems,  the  channel  is  modeled 
and  post-detection  processing  is  developed  to  conform  with 
the  channel  model.  Given  the  prior  probabilities  and  the 
count  parameters  for  Mq  and  ,  and  the  quantization  level, 
J,  the  probability  of  error  for  the  system  can  be  calculated 
directly.  Figure  7  shows  a  simple  flowchart  used  to  calcu¬ 
late  probability  of  error.  Appendix  A,  Figure  15  is  a  list¬ 
ing  of  the  BASIC  language  computer  program  used  to  compute 
data.  Appendix  B,  Table  II  lists  calculated  data  for  both 
quantized  and  unquantized  output,  ideal  direct  detection 
systems . 


INPUTS:  Prior  1 

Number 
Averag* 
Averag< 

3robabilities 
of  Quantizations 
?  Count  for  Mq 
s  Count  for 

_ : _ s 

i _ 

Calculate  N{ 

j  — >  Eq  (29) 

_ 5 

! _ 

Calculate  ( 1— i 

U-I 

5Q)  -**Eq  (30) 
ix) — > Eq  (31) 

_ 2 

> 

Calculate  Lj 

j— >Eq  (36) 

_ ' 

r 

Calculate  P(E)— ->Eq  (37) 

OUTPUTS : 


”0 

<1“P0> 

L0 

P(E) 


d-Pi) 


Figure  7.  Ideal  Detector,  Probability  of  Error  Flowchart 
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IV  Cutoff  Rate  Analysis 


Introduction 

The  design  criterion  for  modulation  systems  in  almost 
universal  use  is  probability  of  error.  While  probability  of 
error  is  meaningful  fdr  "hard  decision"  systems,  where  the 
number  of  outputs  correspond  exactly  to  the  number  of 
inputs,  the  meaning  of  error  in  a  "soft  decision"  system, 
where  the  number  of  outputs  is  greater  than  the  number  of 
inputs,  is  rather  nebulous.  For  this  reason,  another  more 
meaningful  design  criterion  would  prove  useful  for  "soft 
decision"  communication  systems. 

Cutoff  Rate  Development 

Wozencraft  and  Kennedy  (Ref  25)  were  the  first  to  sug¬ 
gest  the  use  of  the  "cutoff  rate,"  Rg,  as  a  design  criterion 
for  a  discrete  memoryless  channel  (DMC) .  The  cutoff  rate 
was  chosen  because  it  is  the  upper  limit  of  code  rates  for 
which  the  average  decoding  computation  per  source  digit  is 
finite  when  sequential  decoding  is  used  (Ref  17:328). 

Massey  (Ref  11)  observed  that  Rg  gives  more  information 
about  the  DMC  than  the  channel  capacity  because  Rg  relates 
both  a  range  of  rates  for  reliable  communication,  as  well  as 
the  necessary  coding  complexity  to  achieve  a  specific  error 
probability. 


As  developed  previously,  a  DMC  of  M  inputs  and  Q 
outputs  can  be  completely  described  by  the  transition 


probabilities,  Pr{q|m},  the  probability  the  channel  output 
is  q  given  the  channel  input  is  m.  The  cutoff  rate,  Rg,  is 
defined  by  using  these  transition  probabilities  as  follows: 


min  Q-l 
P  l 
q=0 


M-l 

Z  p(m) 
m=0 


(39) 


where  p  is  a  probability  distribution  for  the  M  channel 
input  signals.  When  p  is  the  uniform  distribution,  the  cut¬ 
off  rate  defined  by  Eq  (39)  is  called  the  symmetric  cutoff 
rate,  Rg,  given  by: 


(40) 


Additionally,  RQ  -  RQ  for  the  binary  input  case  where  the 
uniform  distribution  is  always  the  minimizing  distribution 
and  in  most  other  cases  of  practical  interest  where  the 
modulation  signal  set  and  the  demodulator  decision  regions 
are  reasonably  symmetric  (Ref  10:438). 


For  convolutional  codes,  Viterbi  (Ref  22)  has  shown  that  the 
decoding  error  probability  is  upperbound  as  follows: 
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(42) 


P ( E )  <  crL2~NR0  if  R  <  Rq 

where 

N  =  code  constraint  length 

L  »  number  of  bits  encoded 

cR  =  multiplying  factor  independent  of  N  and  L 

1 

Rq  =  cutoff  rate 
Throughout  this  report,  block  codes  have  been  used,  thus  Eq 
(41)  applies.  As  stated  previously,  Rq  =  Rq  for  the 
binary  input  case,  so  Eq  (41)  can  be  rewritten: 

P ( E )  <  2“N(R0~R)  if  R  <  Rq  (43) 

Examining  Eq  (43),  it  should  be  noted  that  the  prob¬ 
ability  of  error  bound  will  become  smaller  as  the  symmetric 
cutoff  rate  increases.  Additionally,  the  bound  may  be  made 
arbitrarily  small  by  increasing  N,  as  long  as  the  R  <  RQ 
condition  is  satisfied. 

Thus,  a  new  decision  criteria,  based  on  cutoff  rates 
for  either  block  or  convolutional  codes  can  be  developed. 

To  decrease  the  probability  of  error  bounds  given  by  Eqs 
(41),  (42)  and  (43),  the  cutoff  rate  should  be  maximized. 


Summary 


For  systems  that  employ  "hard  decision"  channels  such 
as  the  unquantized  output  of  an  ideal  unity  gain  detector 
discussed  in  Chapter  II,  the  probability  of  error  design 
criterion  is  quite  meaningful.  However,  if  a  "soft 
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decision"  channel  is  used,  such  as  in  the  quantized  output 
systems  discussed  in  Chapter  III,  the  probability  of  error 
is  meaningless,  and  another  criterion  must  be  developed. 

The  cutoff  rate  criterion  has  proven  to  be  a  very  effective 
method  for  design  of  these  "soft  decision"  channel  systems. 

The  symmetric  cutoff  rate,  Rq,  is  calculated  using  the 
transition  probabilities  that  characterize  the  channel 
model.  This  value  for  Rq  can  then  be  used  to  determine  an 
upper  bound  on  the  probability  of  error  for  the  DMC  given  a 
specific  code  rate  R. 

Appendix  A,  Figure  17  is  a  listing  of  the  BASIC  lan¬ 
guage  computer  program  used  to  calculate  the  symmetric 
cutoff  rate,  Rq,  for  the  various  systems  from  Appendix  B, 

«N» 

Tables  II  and  IV.  Calculated  values  of  RQ  are  presented  in 
Appendix  B,  Tables  III  and  V,  along  with  calculated  P(E) 
values  for  comparison. 
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V 


Avalanche  Detector  with  Output  Quantization 

Introduction 

Chapter  III  examined  optical  direct  detection  receivers 
assuming  an  ideal  gain  factor  of  one.  This  chapter  assumes 
the  receiver  now  uses1 an  avalanche  detector.  Since  the  gain 
of  an  avalanche  detector  is  random  (Ref  15:161),  its  statis¬ 
tics  will  have  an  impact  on  system  performance.  To  determine 
the  impact  that  the  random  gain  has,  the  probability  of  the 
number  of  output  photoelectrons,  given  the  number  of  photo¬ 
electrons  input  to  the  avalanche  detector  is  required. 

While  several  authors  have  conceived  differing  values  for 
the  probability  distribution  of  interest  (Refs  4,14,16),  the 
most  complete  development  has  been  formulated  by  McIntyre 
(Ref  12).  The  probability  distribution  derived  by  McIntyre 
will  be  used  in  this  section  to  provide  the  required  infor¬ 
mation  about  the  avalanche  detector  random  gain.  The 
probability  is  given  by: 


(44) 

where 

n  »  number  of  counts  input  to  the  avalanche  detector 
n+r  «  number  of  counts  output  by  the  avalanche  detector 


Detector 


Figure  8.  Avalanche  Detector  System  Block  Diagram 


k  ■  ratio  of  the  hole  ionization  coefficient  to  the 

electron  ionization  coefficient  for  the  avalanche 
detector 

X  *  n  +  r  -  nM 
nM 

and  M  is  the  average  gain  of  the  avalanche  detector. 

By  using  the  probability  given  by  Eq  (44)  to  modify  the 
equations  from  Chapter  III,  the  effect  of  introducing  an 
avalanche  detector  into  the  optical  communication  system 
becomes  a  logical  extension  of  the  ideal  unity  gain  detector 
previously  examined. 

Channel  Characterization 

The  procedure  used  to  model  the  channel  for  a  communi¬ 
cation  system  with  an  avalanche  detector  parallels  closely 
to  that  used  previously  in  Chapter  III  where  unity  gain  was 
assumed.  Figure  8  shows  a  block  diagram  of  an  avalanche 
detector  system.  Instead  of  each  photon  having  a  count  of 


one,  a  random  value,  g^,  must  be  multiplied  by  each  count  to 
determine  the  total  count  for  time  interval  T.  That  is, 
while  for  the  ideal  unity  gain  detector,  the  total  count  was 
given  by: 

N(T)  =  l  n.  (45) 

i  1 

the  inclusion  of  the  avalanche  detector  results  in  a  total 
count : 

N_  ( T )  =  l  g.n.  (46) 

a  ^  11 

where  g^  is  the  random  gain  associated  with  the  ith  count. 
Thus,  to  make  a  decision  based  on  N  (T),  some  knowledge  of 
the  random  gain,  g^,  should  be  known,  specifically,  the  gain 
probability  distribution  given  by  Eg  (44). 

Repeating  the  logic  used  to  determine  the  effect  of 
output  quantization  developed  in  Chapter  III,  the  count 
N  (T)  can  be  thought  of  as  the  sum  of  J  individual  count 
terms,  N  . .  Again,  to  characterize  the  channel,  the  transi- 
tion  probabilities  for  all  of  the  possible  channel  output 
strings  are  required.  The  output  from  the  channel  will 
again  be  one  of  the  2J  possible  J  bit  long  strings,  where  J 
is  the  number  of  segments  that  T  is  partitioned  into.  Again, 
the  unquantized  system  occurs  when  J  equals  one. 

For  each  of  the  J  segments  of  T,  the  HAP  decision  rule 
of  Eq  (26)  will  be  used  to  determine  whether  the  bit  for 
that  segment  is  a  zero  or  a  one.  Substituting  N  .  for  N.  in 

ai,  1 

Eq  (26)  gives: 
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Returning  to  the  channel  model,  substituting  Eqs  (48) 
and  (49)  into  Eq  (47)  gives: 


n-1 

l  p(t|M.)P  Dx 

t«l  1  >A 


Pr{M0) 


tlx  P^IVPt,n  D0 


TrTM p" 


(53) 


The  probability  of  t  counts  given  with  T  portioned  into  J 
segments  is  given  by  Eq  (27).  Similarly,  the  probability  of 
t  counts  given  Mg  and  J  is  given  by  Eq  (28).  Now,  while 
Eq  (53)  is  not  easily  solved  directly,  it  can  be  solved  by 
iteration,  yielding  a  value,  NQa.  This  NQa  is  analogous  to 
Ng  developed  in  Chapter  III  for  the  ideal  unity  gain  detec¬ 
tor.  Thus,  if  the  count  for  any  segment,  n^,  is  greater 
than  Nga,  that  segment's  bit  is  set  to  one;  if  n^  is  less 
than  Nga ,  the  bit  for  the  segment  is  set  to  zero;  and  again, 
Eq  (18)  is  used  if  n^  *  N0a*  As  An  c^aPter  111 •  this  com¬ 
parison  is  made  for  each  of  the  J  segment  counts,  n^,  until 
a  J  bit  long  string  of  zeros  and  ones  is  formulated.  This 
string  is  the  output  of  the  channel  model,  and,  once  again, 
a  "soft  decision”  is  made. 

To  complete  the  characterization  of  the  channel,  the 
single  bit  crossover  probability  must  be  determined  to 
enable  the  calculation  of  the  final  transition  probabilities 
These  bit  crossover  probabilities  for  each  segment  are  given 


Pr{D1  JMq}  *  Pr{n.  >  NgjMg} 


«  n- 1 

Z 


N 


tJ0  P<tl"o>pt.n  ’  a  ‘  p05 


(54) 


and 


Pr{D0|M1>  =  Pr(ni  <  NQ  jM^ 

SL 


No  -1 


N, 


Z  sa  p(t|M  )P  =  (1_-  p.) 
t-0  n=t+l  1  z*n  1 


(55) 


probabilities  given  by  Eqs  (54)  and  (55)  are  analogous 
±o  the  probabilities  given  by  Eqs  (30)  and  (31)  for  the 
— ^i&deal  unity  gain  detector  system.  -  - 

Por  an  arbitrary  J,  the  transition  probability  associ¬ 
ated  with  a  particular  J  bit  long  string- is  calculated  by 
“*he  procedure  given  in  Chapter  III: 

1  -  Let  the  number  of  ones  in  the  J  bit  string  equal  v, 

then  the  number  of  zeros  in  the  string  is  J  -  v; 

2  -  To  determine  the  transition  probability  of  the 

particular  string,  given  was  transmitted: 


Pr((J  bit  string)^}  ■  P1V(1-P1)J_V  (32) 


where  ( 1— >  is  given  by  Eq  (55); 

3  -  To  determine  the  transition  probability  of  the 
particular  string,  given  MQ  was  transmitted: 
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A  model  which  fully  characterizes  the  channel  employing 
an  avalanche  detector  has  now  been  developed.  The  next  step 
is  to  determine  the  optimum  post-detection  processing 
required  for  the  channel  model. 


Post-Detection  Processing 

The  MAP  decision  rule  of  Eq  (3)  is  once  again  the  basis 
for  developing  the  optimum  post-detection  processor.  For 
any  J  bit  string  output  by  the  channel  the  MAP  decision  rule 
Ogives:"-  —  -  -------  ■ 

jl  Pr{M0) 

<  PrTh.} 

D0  (34) 

where  X^X2..«XJ  is  the  J  bit  long  string  of  zeros  and  ones. 
Substituting  Eqs  (32)  and  (33)  into  Eq  (34)  givec: 


■ao-E  -Tcr-rv 


Pr{X1X2, 


TFOcTX 


.XjjMj} 


1~2 


•W 


iTZ-Zl  - 


v,.  „  -  v  D. 

(1  ‘  Pj) 

”rP0  ^  '  Po^  d, 


Pr{M0) 

PFTHJT 


(35) 


where  v  is  the  number  of  ones  in  the  J  bit  string.  This  is 
identical- to  the- results  developed  in  Chapter  III  for  the 
ideal  unity  gain  detector.  Similarly,  Ln  as  defined  by: 


1 


In 


will  be  used  as  the  decision  maker  of  Figure  6.  Thus,  the 

t 

major  difference  between  the  avalanche  detector  system  and 
the  ideal  unity  gain  detector  system  is  in  the  calculation 
of  the  transition  probability  for  each  of  the  J  individual 
bits,  as  is  evident  by  comparing  Eq  (54)  with  Eq  (30)  and 
Bq  (55)  with  Eq  (31). 


Probability  of  Error 

The  system  probability  of  error  is  given  by: 


P(E)  -  Pr{D0|M1)Pr{M1}  +  Pr{D1|M0)Pt{M0>  (5) 

which  reduces  to: 


P(E) 


M  /j\  T 

MvK  (1 ' Pl)' 


Pr{Mj} 


v 


(1  -  p0)v 


Pr{M0) 


(37) 


as  was  the  case  for  the  ideal  unity  gain  detector. 

Summary 

The  addition  of  an  avalanche  detector  into  the  system 
adds  a  random  factor  into  the  computations  and  thus  adds 


complexity  to  the  channel  model.  By  using  the  probability 

distribution  for  PR  n+r  developed  by  McIntyre,  the  channel 

can  be  characterized,  and  probability  of  error  data  is 

readily  calculated.  [It  should  be  noted  that  a  discrepancy 

was  noted  in  the  paper  by  McIntyre  (Ref  12).  Using  Eq  (44) 

to  calculate  P  ,  ,  the  results  obtained  did  not  agree  with 
n,n+r 

Figure  3  of  McIntyre's  work.  If,  however,  the  results  of 

Eq  (44)  were  multiplied  by  n.  Figure  3  was  reproduced.  For 

the  purposes  of  this  report,  the  results  of  Eq  (44)  were 

used  in  the  calculations  where  P„  were  required.] 

n,n+r 

Figure  9  shows  a  simple  flowchart  used  in  calculating 
the  system  probability  of  error.  Appendix  A,  Figure  19  is  a 
listing  of  the  BASIC  language  computer  program  used  to  com¬ 
pute  data.  Appendix  B,  Table  IV  lists  calculated  data  for 
the  avalanche  detector  system. 


INPUTS:  Prior  Probabilities 

Number  of  Quantizations 
Average  Counts  for  and  M, 
Average  Gain  of  Avalanche  Detector 
Ionization  Coefficient  Ratio 


Calculate  (l-pQ) — >Eq  (54) 
(1-p,) — »Eq  (55) 


Calculate  LQ — ^Eq  (36) 


Calculate  P(E) — >Eq  (37) 


OUTPUTS: 


U-P0>  (l-P^ 


P(E) 


Figure  9.  Avalanche  Detector,  Probability  of  Error 

Flowchart 


VI  Conclusions 


General  Information 

Several  equations  have  been  developed  to  assist  in 
analyzing  the  Multilevel  Amplitude  Keyed  direct  detection 
optical  communication' system.  The  discrete  nature  of  the 
integrate-and-dump  receiver  and  the  counting  process  and  how 
the  equations  are  affected  by  the  discreteness  lead  to  some 
important  conclusions. 

To  characterize  the  ideal  unity  gain  detector  channel, 
the  threshold  Nn  _  was  defined  to  be: 


(29) 


This  threshold  determined  whether  the  bit  for  a  segment  was 
to  be  a  zero  or  a  one.  It  is  important  to  note  that  for 
Ng  j  <  0  ,  the  bit  will  always  be  set  to  one;  M1  will  always 
be  decided;  and  the  system  probability  of  error  will  be  the 
prior  probability  of  Mq.  This  will  occur  for  MQ  <  M^  , 

Ug  <  ,  and  large  J.  For  Mg  >  M^  and  ,  this 

situation  is  reversed;  Nn  _  is  positive  and  more  bits  will 

U  f  J 

be  set  to  zero.  Eventually,  for  mq>>m1  ,  and  small  J,  all 
the  bits  will  be  set  to  zero;  Mg  will  always  be  decided;  and 
the  system  probability  of  error  will  equal  the  prior  prob¬ 
ability  of  M. .  The  second  important  fact  concerning  Nn  .  is 


that  a  threshold  of  2.99  will  give  equivalent  results  as  a 
threshold  of  2.01.  Clearly,  due  to  the  discrete  nature  of 
the  counting  process,  some  information  about  the  relative 
position  between  a  threshold  and  the  count  is  lost. 

A  similar  situation  exists  for  Lq,  the  threshold  for 
both  ideal  unity  gain* and  avalanche  detector  systems,  given 
by: 


(36) 

If  Lq  £  0  ,  M1  is  always  decided,  and  the  probability  of 
error  equals  the  prior  probability  of  MQ.  If  LQ  >  J  ,  MQ 
is  always  decided,  and  the  probability  of  error  equals  the 
prior  probability  of  .  Likewise,  a  threshold  of  2.99  will 
result  in  the  same  decision  as  a  threshold  of  2. 011  Thus, 

Nq  j  and  Lq  are  both  affected  in  a  similar  manner  by  the 
discrete  nature  of  the  counting  process. 

The  final  equation  of  interest  is: 


p(tlMl)Pt,n  °1  Pr (Mq) 

n-1  <  Pr{M. } 

:0>Pt,n  D0 


(53) 


This  equation  can  be  solved  by  iteration  for  n,  given  a 
value  Nqs.  The  NQa  is  a  threshold,  comparable  to  NQ  j,  but 
is  used  in  analysis  of  an  avalanche  detector  system.  This 


equation  may  not  yield  a  threshold  for  certain  values  of 
prior  probabilities  and  average  counts  associated  with  MQ 
and  M^.  For  example,  given  equilikely  priors,  Eq  (53) 
reduces  to: 


t;;  p^i^t.n  “i 


Since  P.  „  will  be  the  same  value  in  both  the  numerator  and 
the  denominator  for  any  value  of  t  and  n,  in  order  to  solve 
Eq  (54)  for  NQa,  p(tjMQ)  must  be  greater  than  p(t|M1)  for 
some  integer  value  of  t.  Using  Eqs  (8)  and  (9),  it  can  be 
shown  that  this  is  not  necessarily  the  case.  [For  example: 
if  ■  .2  and  y1  »  1  ,  then  p(t|MQ)  <  p(t|M^)  for  all 
integer  values  of  t.]  Thus,  care  must  be  used  in  selecting 
values  for  yQ  and  y1  when  analyzing  the  avalanche  detector 
system. 


Ideal  Detector,  Quantized  Output 

Using  the  calculated  data  from  Table  II,  Appendix  B, 
two  distinct  conclusions  concerning  probability  of  error  can 
be  made.  First,  as  shown  in  Figure  10,  given  yQ,  the 
average  count  for  MQ,  as  y^^  gets  further  from  yQ,  the  prob¬ 
ability  of  error  decreases.  This  is  as  expected  because  as 
the  "distance"  between  two  signals  increases,  the  probabil¬ 
ity  of  making  an  error  should  decrease.  The  second 
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Figure  10.  Probability  of  Error  vs  Quantizations  for  Ideal  Detector 

(yQ  constant  :  y,  varying) 


] 


Figure  11.  Probability  of  Error  vs  Quantizations  for  Ideal  Detector 

(yn,  y,  varying) 


conclusion  is  not  as  obvious.  As  shown  in  Figure  11,  as  the 
quantizations  increase,  the  probability  of  error  first 
increases  and  then  decreases  asymptotically  to  the  level  for 
unquantized  output.  The  reason  for  this  is  that  by  quantiz¬ 
ing  the  output,  a  "hard  decision"  channel  became  a  "soft 
decision"  channel  where  probability  of  error  is  no  longer  a 
valid  decision  criterion. 

Cutoff  Rate  Analysis 

The  conclusion  made  by  examining  the  cutoff  rates  has 
already  been  stated  by  examining  Eq  (41).  As  shown  in 
Figure  12,  for  the  ideal  unity  gain  detector,  and  Figure  13, 
for  the  avalanche  detector,  as  the  symmetric  cutoff  rate  RQ 
increases,  the  probability  of  error  decreases.  Thus,  while 
probability  of  error  has  no  meaning  for  "soft  decision" 
channels,  examination  of  the  cutoff  rate,  specifically  maxi¬ 
mizing  the  cutoff  rate,  will  result  in  a  criterion  that  can 
effectively  be  used  to  minimize  probability  of  error. 

Avalanche  Detector,  Quantized  Output 

Plotting  the  data  from  Table  IV,  Appendix  B,  in  Figure 
14  shows  that  the  minimum  probability  of  error  is  not 
obtained  with  unquantized  output,  as  it  was  for  the  ideal 
unity  gain  detector.  The  value  of  J  that  returns  minimum 
probability  of  error  increases  as  the  distance  between  Uq 
and  increases.  Additionally,  the  value  of  J  that  pro¬ 
duces  minimum  probability  of  error  also  results  in  Npa  being 
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Figure  12.  Probability  of  Error  vs  Symmetric  Cutoff  Rate  for  Ideal  Detector 


Figure  13.  Probability  of  Error  vs  Symmetric  Cutoff  Rate  for  Avalanche  Detector 


J  (Quantization  Levels) 

Figure  14.  Probability  of  Error  vs  Quantizations  for  Avalanche  Detectors 


a  minimum  value  and  LQ  being  a  maximum  value.  The  data  from 
Table  IV  is  in  agreement  with  the  ideal  unity  gain  detector 
data  in  that  as  the  distance  between  Uq  and  increases, 
the  probability  of  error  decreases,  as  expected. 

Topics  for  Future  Research 

i 

An  obvious  extension  of  this  paper  would  be  to  abandon 
the  MAP  decision  rule  and  try  to  improve  the  probability  of 
error  results  by  using  a  cutoff  rate  criterion  when  the  out¬ 
put  is  quantized  and  a  "soft  decision"  channel  exists.  The 
cutoff  rate  criterion  results  could  be  compared  to  the 
results  obtained  by  the  MAP  criterion  of  this  paper  to 
demonstrate  the  benefits  of  using  cutoff  rate  criterion  for 
•soft  decisions." 

Another  area  of  interest  would  be  expanding  the  study 
of  the  avalanche  detector  system  and  attempting  to  determine 
what  level  of  quantization  would  return  a  minimum  probabil¬ 
ity  of  error  for  any  given  MQ  and  M1  characteristics. 

These  are  just  two  areas  barely  touched  in  this  work 
that  may  provide  additional  depth  of  knowledge  into  the 
field  of  optical  communications. 
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Figure  15  is  the  listing  of  the  BASIC  computer  program 
used  to  calculate  the  data  listed  in  Table  II  of  Appendix  B. 

A  sample  of  the  output  from  this  program  is  given  in  Figure 
16.  The  program  requires  the  prior  probabilities  of  messages 
Mq  and  the  number  of  segments  in  the  interval  T  and  the 
average  count  for  messages  Mg  and  as  input  variables.  As 
output#  the  program  returns  NQ  [calculated  by  Eq  (29)],  the 
bit  transition  probabilities,  [calculated  by  Eqs  (30)  and 
(31)],  Lg  [calculated  by  Eq  (36)],  and  the  system  probability 
of  error  [calculated  by  Eq  (37)]. 
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LPRINT 

LPRINT 


LPRINT 

LPRINT 

LPRINT 

LPRINT 

LPRINT 


LPRINT 

LPRINT 

LPRINT 

LPRINT 

LPRINT 

LPRINT 

LPRINT 

LPRINT 


"THESIS1  -  THIS  PROGRAM  CALCULATES:" 

"  N0  -  THE  NUMBER  OF  COUNTS  PER  SEGMENT  USED 

TO  DETERMINE  A  1  OR  0  FOR  THAT 
SEGMENT" 

"  IF  N  >  N0  DECIDE  1  FOR  THAT  SEGMENT" 

"  IF  N  <  N0  DECIDE  0  FOR  THAT  SEGMENT" 

"  CPI  -  THE  PROBABILITY  OF  DECIDING  A  0  FOR  A 

SEGMENT  GIVEN  MESSAGE  1  WAS  SENT" 

"  CP0  -  THE  PROBABILITY  OF  DECIDING  A  1  FOR  A 

SEGMENT  GIVEN  MESSAGE  0  WAS  SENT" 

"  L0  -  THE  NUMBER  OF  SEGMENTS  WHERE  A  1  WAS 

DECIDED  REQUIRED  TO  MAKE  THE  FINAL 
DECISION  THAT  MESSAGE  1  WAS  SENT" 

"  IF  MORE  THAN  L0  OF  THE  J  SEGMENTS  ARE 

1  THEN  DECIDE  MESSAGE  1  WAS  SENT" 

"  IF  LESS  THAN  L0  OF  THE  J  SEGMENTS  ARE 

1  THEN  DECIDE  MESSAGE  0  WAS  SENT" 

"  P  OF  E  -  THE  OVERALL  SYSTEM  PROBABILITY  OF 

ERROR 

"GIVEN  THE  INPUTS:" 

"  PR0  -  THE  PRIOR  PROBABILITY  OF  MESSAGE  0" 

"  J  -  THE  NUMBER  OF  SEGMENTS  T  IS  DIVIDED 

INTO" 

"  M0  -  THE  AVERAGE  COUNT  FOR  MESSAGE  0  BEING 
SENT’ 

"  Ml  -  THE  AVERAGE  COUNT  FOR  MESSAGE  1  BEING 
SENT" 


LPRINT 

LPRINT 

LPRINT 

LPRINT  "  J  PR0  PR1  M0  Ml",,"  N0","  CP0" , " 

CPI’,"  L0" , "  P  OF  E" 

INPUT  "ENTER  PRIOR  PROBABILITY  OF  MESSAGE  0";  PR0 
INPUT  "ENTER  NUMBER  OF  SEGMENTS  IN  T";  J 
INPUT  "ENTER  AVERAGE  COUNT  FOR  MESSAGE  0";  M0 
INPUT  "ENTER  AVERAGE  COUNT  FOR  MESSAGE  1";  Ml 
PR1  -  1  -  PR0 
REM 

REM  CALCULATE  N0 

REM 

N0  =  ((LOG  ( PR0/PR1 ) )  +  (M1-M0)/J)  /  LOG(M1/M0) 

IF  N0<=0  GOTO  780 
REM 

REM  CALCULATE  CP0  AND  CPI 

REM 

SM0  =  0 
SMI  =  0 

FOR  N  =  0  TO  INT(N0) 


Figure  15.  Ideal  Detector,  Quantized  Output  Program  Listing 
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TERM#  =  (  M0 / J  )  a  N 
TER  Ml  =  (M1/J)AN 
FOR  A  »  1  TO  N 
TERM#  =  TERM#/ A 
TERMl  =  TERM1/A 
NEXT  A 

SM#  =  SM#  +  TERM# 

SMI  =  SMI  +  TERMl 
NEXT 

CP#  =  1  -  (SM#*(EXP(-M#/J) ) ) 

CPI  =  SMI* ( EXP { -Ml/J } ) 

REM 

REM  CALCULATE  L# 

REM 

L#  =  ( (LOG(PR#/PRl) )-  J* ( LOG ( CPI/ ( 1-CP0) ) ) ) / 

LOG ( ( (1-CP0)*{1-CP1) )/(CP#*CPl) ) 

REM 

REM  CALCULATE  P  OF  E 

REM 

SUM1  =  0 

FOR  R  =  0  TO  INT { L# ) 

GOSUB  67# 

SUM1  =  SUM1  +  ( NUMBER ) * ( CPI a( J-R ) ) * ( 1-CPl )aR 
NEXT 

SUM#  =  # 

FOR  R  =  INT ( L0+1 )  TO  J 
GOSUB  67# 

SUM#  =  SUM#  +  (NUMBER) *( (1-CP#)a( J-R) )*CP#aR 
NEXT 

ERRPROB  =  SUM1 *PRl  +  SUM#*PR0 

LPRINT  J ; "  " ; PR# ; "  ";PRl;"  ";M0;"  " ; Ml , , N0, CP#, CPI , 

L#,  ERRPROB 

END 

NUMBER  =  1 

IF  J  =  R  GOTO  77# 

B  =  J 
C  =  J  -  R 

NUMBER  =  NUMBER *B/C 
B  =  B  -  1 
C  *  C  -  1 

IF  B  -  R  THEN  B  =  1 
IF  C  ■  #  THEN  C  =  1 
IF  B  >  C  GOTO  710 
RETURN 

IF  PR#<PR1  THEN  CP0=1 : CP1=# : L0=# : ERRPROB=PR#  ELSE 
CP0=0 : CP1=1 : L#=J : ERRPROB=PRl 
GOTO  65# 
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THESIS1  -  THIS  PROGRAM  CALCULATES: 

N0  -  THE  NUMBER  OF  COUNTS  PER  SEGMENT  USED  TO  DETERMINE 
A  1  OR  0  FOR  THAT  SEGMENT 
IF  N  >  N0  DECIDE  1  FOR  THAT  SEGMENT 

IF  N  <  N0  DECIDE  0  FOR  THAT  SEGMENT 

CPI  -  THE  PROBABILITY  OF  DECIDING  A  0  FOR  A  SEGMENT 
GIVEN  MESSAGE  1  VJAS  SENT 

CP0  -  THE  PROBABILITY  OF  DECIDING  A  1  FOR  A  SEGMENT 
GIVEN  MESSAGE  0  WAS  SENT 

L0  -  THE  NUMBER  OF  SEGMENTS  WHERE  A  1  WAS  DECIDED 

REQUIRED  TO*  MAKE  THE  FINAL  DECISION  THAT  MESSAGE  1 
WAS  SENT 

IF  MORE  THAN  L0  OF  THE  J  SEGMENTS  ARE  1  THEN 
DECIDE  MESSAGE  1  WAS  SENT 

IF  LESS  THAN  L0  OF  THE  J  SEGMENTS  ARE  1  THEN 
DECIDE  MESSAGE  0  WAS  SENT 
P  OF  E  -  THE  OVERALL  SYSTEM  PROBABILITY  OF  ERROR 
GIVEN  THE  INPUTS: 

PR0  -  THE  PRIOR  PROBABILITY  OF  MESSAGE  0 

J  -  THE  NUMBER  OF  SEGMENTS  T  IS  DIVIDED  INTO 


M0 

Ml 

-  THE 

-  THE 

AVERAGE  COUNT 
AVERAGE  COUNT 

FOR  MESSAGE 
FOR  MESSAGE 

0  BEING 
1  BEING 

SENT 

SENT 

J 

PR0 

PR1 

M0 

HI 

N0 

10 

.5 

.5 

5 

25 

1.24267 

CP0 

CPI 

L0 

P  OF  E 

.090204 

.287297 

3.58016 

8.54147E-03 

Figure  16.  Ideal  Detector,  Quantized  Output  Program 

Sample  Output 


Cutoff  Rate  Calculations 

Figure  17  is  the  listing  of  the  BASIC  computer  program 
used  to  calculate  the  data  listed  in  Table  III  of  Appendix  B. 
A  sample  of  the  output  from  this  program  is  given  in  Figure 
18.  The  program  requires  the  bit  transition  probabilities 
(calculated  by  Eqs  (27)  and  (28)]  and  the  number  of  segments 
in  the  interval  T  as  input  variables.  The  program  returns 
the  symmetric  cutoff  rate  [calculated  by  Eq  (40)]  and  the 
system  probability  of  error  bound  [calculated  by  Eq  (43)]  as 
output. 
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10 

LPRINT  "THESIS2  -  THIS  PROGRAM  CALCULATES:" 

20 

LPRINT  "  R0  -  THE  SYMMETRIC  CUT-OFF  RATE" 

40 

LPRINT  "GIVEN  THE  INPUTS:" 

50 

LPRINT  "  CPI  -  THE  PROBABILITY  OF  DECIDING 

A  0  FOR 

A 

SEGMENT  GIVEN  MESSAGE  1  WAS 

SENT" 

60 

LPRINT  "  CP0  -  THE  PROBABILITY  OF  DECIDING 

A  1  FOR 

A 

SEGMENT  GIVEN  MESSAGE  0  WAS 

SENT" 

70 

LPRINT  "  J  -  THE  NUMBER  OF  SEGMENTS  T  IS 

DIVIDED 

INTO" 

80 

LPRINT 

90 

LPRINT 

100 

LPRINT 

110 

LPRINT  "  J  CP0  CPI","  R0" 

120 

INPUT  "ENTER  NUMBER  OF  SEGMENTS  IN  T";  J 

130 

INPUT  "ENTER  CROSSOVER  PROBABILITY  FOR  H0";  CP0 

140 

INPUT  "ENTER  CROSSOVER  PROBABILITY  FOR  HI";  CPI 

150 

RSUM  =  0 

160 

REM 

170 

REM  CALCULATE  R0 

180 

REM 

190 

FOR  L  =  0  TO  J 

200 

GOSUB  310 

210 

MSUM  •  NUMBER  *  ( SQR ( ( 1-CP0 ) a( J-L )  *(CP0aL))  + 

SQR ( ( CP1A (J-L) )*(1-CP1)AL) )A  2 

220 

RSUM  =  RSUM  +  MSUM 

230 

NEXT 

240 

R0  -  1  -  ( (LOG(RSUM/2))/LOG  (2)) 

290 

LPRINT  J;"  " ;CP0; "  ";CP1;"  ";R0 

300 

END 

310 

NUMBER  =  1 

320 

IF  J  =  L  GOTO  410 

330 

B  =  J 

340 

C  =  J  -  L 

350 

NUMBER  =  NUMBER*B/C 

360 

B  =  B  -  1 

370 

C  =  C  -  1 

380 

IF  B  =  L  THEN  B  =  1 

390 

IF  C  =  0  THEN  C  =  1 

400 

IF  B  >  C  GOTO  350 

410 

RETURN 

Figure  17.  Cutoff  Rate  Calculations  Program  Listing 
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THESIS2  -  THIS  PROGRAM  CALCULATES: 

R0  -  THE  SYMMETRIC  CUT-OFF  RATE 
GIVEN  THE  INPUTS: 

CPI  -  THE  PROBABILITY  OF  DECIDING  A  0  FOR  A  SEGMENT 
GIVEN  MESSAGE  1  WAS  SENT 

CP0  -  THE  PROBABILITY  OF  DECIDING  A  1  FOR  A  SEGMENT 
GIVEN  MESSAGE  0  WAS  SENT 
J  -  THE  NUMBER  OF  SEGMENTS  T  IS  DIVIDED  INTO 


CP0  . 


,32968 


.135335 


.508227 


Figure  18.  Cutoff  Rate  Calculations  Program  Sample  Output 


Avalanche  Detector ,  Quantized  Output 

Figure  19  is  the  listing  of  the  BASIC  computer  program 
used  to  calculate  the  data  listed  in  Table  IV  of  Appendix  B. 

A  sample  of  the  output  from  this  program  is  given  in  Figure 
20.  The  program  requires  the  prior  probabilities  of  messages 
Mg  and  M^,  the  number- of  segments'  in  the  interval  T,  the 
average  count  for  messages  Mg  and  M^,  the  average  gain  of 
the  avalanche  detector  and  the  ionization  coefficient  ratio 
of  the  avalanche  detector  as  input  variables.  The  output 
consists  of  Nga  [calculated  from  Eq  (53)],  the  bit  transition 
probabilities  [calculated  from  Eqs  (54)  and  (55)],  Lg 
[calculated  by  Eq  (36)],  and  the  system  probability  of  error 
[calculated  by  Eq  (37)]. 


10 

LPRINT 

20 

LPRINT 

30 

LPRINT 

40 

LPRINT 

50 

LPRINT 

60 

LPRINT 

70 

LPRINT 

80 

LPRINT 

90 

LPRINT 

100 

LPRINT 

110 

LPRINT 

120 

LPRINT 

130 

LPRINT 

140 

LPRINT 

150 

LPRINT 

160 

LPRINT 

170 

LPRINT 

180 

LPRINT 

190 

LPRINT 

200 

LPRINT 

210 

INPUT 

220 

INPUT 

230 

INPUT 

240 

INPUT 

250 

INPUT 

260 

INPUT 

270 

GOSUB 

280 

REM 

290 

REM 

300 

REM 

310 

SUM0  = 

320 

CPI  = 

"THESIS3  -  THIS  PROGRAM  CALCULATES; 

"  MO  -  THE  NUMBER  OF  COUNTS  PER  SEGMENT  USED 

TO  DETERMINE  A  1  OR  0  FOR  THAT 
SEGMENT" 

"  IF  N  >  MO  DECIDE  1  FOR  THAT  SEGMENT" 

"  IF  N  <  MO  DECIDE  0  FOR  THAT  SEGMENT" 

■  CPI  -  THE  PROBABILITY  OF  DECIDING  A  0  FOR  A 
SEGMENT  GIVEN  MESSAGE  1  WAS  SENT" 

"  CP0  -  THE  PROBABILITY  OF  DECIDING  A  1  FOR  A 
SEGMENT  GIVEN  MESSAGE  0  WAS  SENT" 

"  L0  -  THE  NUMBER  OF  SEGMENTS  WHERE  A  1  WAS 

DECIDED  REQUIRED  TO  MAKE  THE  FINAL 
DECISION  THAT  MESSAGE  1  WAS  SENT" 

"  IF  MORE  THAN  L0  OF  THE  J  SEGMENTS  ARE 

1  THEN  DECIDE  MESSAGE  1  WAS  SENT" 

"  IF  LESS  THAN  L0  OF  THE  J  SEGMENTS  ARE 

1  THEN  DECIDE  MESSAGE  0  WAS  SENT" 

"  P  OF  E  -  THE  OVERALL  SYSTEM  PROBABILITY  OF 
ERROR" 

"GIVEN  THE  INPUTS: " 

"  PR0  -  THE  PRIOR  PROBABILITY  OF  MESSAGE  0" 

"  J  -  THE  NUMBER  OF  SEGMENTS  T  IS  DIVIDED 

INTO" 

"  M0  -  THE  AVERAGE  COUNT  FOR  MESSAGE  0  BEING 

SENT" 

"  Ml  -  THE  AVERAGE  COUNT  FOR  MESSAGE  1  BEING 

SENT" 

"  M  -  THE  AVERAGE  GAIN  OF  THE  AVALANCHE 

DETECTOR" 

"  K  -  THE  HOLE  IONIZATION  COEFFICIENT  TO 

ELECTRON  IONIZATION  COEFFICIENT  RATIO" 


"ENTER  PRIOR  PROBABILITY  OF  MESSAGE  0";  PR0 

"ENTER  NUMBER  OF  SEGMENTS  IN  T";  J 

"ENTER  AVERAGE  COUNT  FOR  MESSAGE  H0";  M0 

"ENTER  AVERAGE  COUNT  FOR  MESSAGE  HI";  Ml 

"ENTER  AVERAGE  GAIN  FOR  THE  AVALANCHE  DETECTOR" ;M 

"ENTER  HOLE  IONIZATION  COEFFICIENT  TO  ELECTRON 

IONIZATION  COEFFICIENT  RATIO";  K 


CALCULATE  CP0  AND  CPI 


Figure  19.  Avalanche  Detector,  Quantized  Output 

Program  Listing 

(Sheet  1  of  3) 
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330  FOR  NI  »  1  TO  INT( MO) -1 
340  FOR  MCT  =  NI+1  TO  INT(MO) 

350  GOSUB  810 

360  TERM0  =  P0(NI)  *  PNMO 

370  TERMl  =  Pl(NI)  *  PNMO 

380  SUM0  =  SUM0  4-  TERM0 

390  CPI  =  CPI  +  TERMl 

400  NEXT  MCT 

410  NEXT  NI 

420  CPI  =  CPI  +  PI (0) 

430  CP0  =  1  -  (SUM0  4  P0 ( 0 ) ) 

440  GOTO  850 
450  REM 

460  REM  CALCULATE  MO 
470  REM 

480  PR1  =  1  -  PR0 

490  DIM  P0(INT(M1*M/J) ) , Pi ( INT( M1*M/J) ) 

500  P0 (0 )  =  EXP ( -M0/ J ) 

510  PI (0)  =  EXP ( -Ml/ J ) 

520  FOR  N  =  1  TO  ( INT(M1*M/J)-1 ) 

530  P0(N)  =  P0(N-1)*(M0/(J*N)) 

540  PI ( N )  =  Pi ( N-l ) * ( Ml/ ( J*N ) ) 

550  NEXT  N 

560  RIGHT  =  PR0/PR1 

570  MO  *  M 

580  DIFFLT  *  1E+37 

590  SUM  =  0 

600  TOPSUM  =  0 

610  BOTSUM  =  0 

620  FOR  NI  ■  1  TO  INT(MO-.001) 

630  MCT  *  MO 
640  GOSUB  810 
650  TOP  =  PI { NI ) *PNMO 
660  BOT  =  P0 ( NI ) *PNMO 
670  TOPSUM  =  TOPSUM  +  TOP 
680  BOTSUM  «  BOTSUM  +  BOT 
690  NEXT  NI 

700  LEFT  *  TOPSUM/BOTSUM 
710  DIFF  *  RIGHT  -  LEFT 

720  IF  ABS(DIFF)  >  ABS (DIFFLT)  GOTO  760  ELSE  DIFFLT  =  DIFF 
730  LMO  =  MO 

740  IF  DIFF  >  0  THEN  MO  =  l.l*MO  ELSE  MO  =  ,9*MO 
750  GOTO  590 

760  MO  »  (DIFFLT/ ( DIFFLT-DIFF ) ) * ( MO-LMO) +LMO 
770  RETURN 
780  REM 

790  REM  CALCULATE  PNMO 


Figure  19,  Avalanche  Detector,  Quantized  Output 
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800  HEM 

810  P  =  MCT/(NI*M) 

820  R  *  MCT  -  NI 

830  PNMO  -  (NI/ ( ( 6 . 28*MCT* ( NI+K*R) *R)A . 5 ) ) * ( ( 1- ( NI* ( F-l/R) ) 
AR)*((1+(NI*(1-K)*(F-1)/(NI+K*R) ) )A ( (NI+K*R)/ 

( 1-K) ) ) 

840  RETURN 
850  REM 

860  REM  CALCULATE  L0 

870  REM 

880  L0  =  ( (LOG(PR0/PR1) )-  J* ( LOG ( CPI/ ( 1-CP0) ) ) )/ 

LOG( ( (1-CP0)*(1-CP1) )/(CP0*CPl) ) 

890  REM 

900  REM  CALCULATE  P  OF  E 

910  REM 

920  SUMl  =  0 

930  FOR  R  =  0  TO  INT(L0) 

940  G0SU3  1060 

950  SUMl  =  SUMl  +  (NUMBER) *(CP1a(J-R) )*(1-CP1)aR 
960  NEXT 
970  SUM0  =  0 

980  FOR  R  =  INT ( L0+1 )  TO  J 
990  GOSUB  1060 

1000  SUM0  =  SUM0  +  ( NUMBER ) * ( ( 1-CP0 ) A ( J-R ) ) *CP0A  R 
1010  NEXT 

1020  ERRPROB  =  SUM1*PR1  +  SUM0*PR0 

1030  LPRINT  "  J  PR0  PR1  M0  Ml  K  M","  MO", 

"  CP0" , "  CPI","  L0" , "  P  OF  E" 

1040  LPRINT  J;"  ";PR0;"  ";PR1;"  ";M0;"  ";M1;"  "? Ky" 

" ;M; ’  ", MO, CP0, CP 1,L0, ERRPROB 

1050  END 

1060  NUMBER  =  1 

1070  IF  J  =  R  GOTO  1160 

1080  B  =  J 

1090  C  =  J  -  R 

1100  NUMBER  =  NUMBER*B/C 

1110  B  =  B  -  1 

1120  C  =  C  -  1 

1130  IF  B  =  R  THEN  B  =  1 

1140  IF  C  =  0  THEN  C  =  1 

1150  IF  B  >  C  GOTO  1100 

1160  RETURN 


Figure  19.  Avalanche  Detector,  Quantized  Output 

Program  Listing 
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THESIS3  - 
MO  - 


CPI  - 


CP0  - 
L0  - 


P  OF  E  - 
GIVEN  THE 
PR0 
J 
M0 
Ml 
M 
K 


THIS  PROGRAM  CALCULATES: 

THE  NUMBER  OF  COUNTS  PER  SEGMENT  USED  TO  DETERMINE 

A  1  OR  0  FOR  THAT  SEGMENT 

IF  N  >  MO  DECIDE  1  FOR  THAT  SEGMENT 

IF  N  <  MO  DECIDE  0  FOR  THAT  SEGMENT 

THE  PROBABILITY  OF  DECIDING  A  0  FOR  A  SEGMENT  GIVEN 

MESSAGE  1  WAS  SENT 

THE  PROBABILITY  OF  DECIDING  A  1  FOR  A  SEGMENT  GIVEN 
MESSAGE  0  WAS  SENT 

THE  NUMBER  OF  SEGMENTS  WHERE  A  1  WAS  DECIDED 
REQUIRED  TO.  MAKE  THE  FINAL  DECISION  THAT  MESSAGE  1 
WAS  SENT 

IF  MORE  THAN  L0  OF  THE  J  SEGMENTS  ARE  1  THEN  DECIDE 
MESSAGE  1  WAS  SENT 

IF  LESS  THAN  L0  OF  THE  J  SEGMENTS  ARE  1  THEN  DECIDE 
MESSAGE  0  WAS  SENT 

THE  OVERALL  SYSTEM  PROBABILITY  OF  ERROR 
INPUTS: 

-  THE  PRIOR  PROBABILITY  OF  MESSAGE  0 

-  THE  NUMBER  OF  SEGMENTS  T  IS  DIVIDED  INTO 

-  THE  AVERAGE  COUNT  FOR  MESSAGE  0  BEING  SENT 

-  THE  AVERAGE  COUNT  FOR  MESSAGE  1  BEING  SENT 

-  THE  AVERAGE  GAIN  OF  THE  AVALANCHE  DETECTOR 

■  THE  HOLE  IONIZATION  COEFFICIENT  TO  ELECTRON 
IONIZATION  COEFFICIENT  RATIO 


P  OF  E 


4.38521 


.475774 


149137 


3.41895 


.160942 


Figure  20.  Avalanche  Detector,  Quantized  Output  Program 

Sample  Output 


Table  II  lists  computed  values  of  Ng,  (1-Pg),  (1-p^), 
Lg,  and  P { E )  for  several  values  of  J,  Ug,  jj  ^  and  equiprob- 
able  priors,  for  the  ideal  unity  gain  detector.  Table  III 
lists  computed  values  of  Rg  given  the  crossover  probabili¬ 
ties,  (1-pg)  and  (1-p^),  from  Table  II.  Table  IV  lists 
computed  values  of  Ng,  (1-Pg),  (1-p^),  Lg,  and  P(E)  for 
several  values  of  J,  Ug,  and  equiprobable  priors,  for 
the  avalanche  detector  system.  An  avalanche  detector  with 
average  gain  of  10  and  ionization  coefficient  ratio  of  .5 
was  used  in  calculating  the  values  of  Table  IV.  These 
values  were  used  only  to  shorten  calculation  time  and  pro¬ 
vide  representative  data.  Table  V  lists  computed  values  of 
Rg  given  the  crossover  probabilities,  (1-pg)  and  (1-Pj), 


from  Table  IV. 


TABLE  II 


Ideal  Detector,  Quantized  Output 
Pr{MQ}  =  Pr{M1> 


J 

^0 

^1 

^0 

U-Pq) 

d-Pi) 

^0 

P(E) 

1 

1 

4 

2.16 

0.080 

0.238 

0.38 

0.1592 

2 

1 

4 

1.08 

0.090 

0.406 

0.60 

0.1686 

3 

1 

4 

0.72 

0.283 

0.264 

1.53 

0.1837 

,4 

1 

4 

0.54 

0.221 

0.368 

1.67 

0.1792 

5 

1 

4 

0.43 

0.181 

0.449 

1.75 

0.1777 

10 

1 

4 

0.22v 

0.095 

0.670 

1.94 

0.1768 

15 

1 

4 

0.14 

0.064 

0.766 

2.01 

0.1750 

20 

1 

4 

0.11 

0.049 

0.819 

2.05 

0.1706 

30 

1 

4 

0.07 

0.033 

0.875 

2.09 

0.1665 

40 

1 

4 

0.05 

0.025 

0.905 

2.11 

0.1645 

50 

1 

4 

0.04 

0.020 

0.923 

2.12 

0.1634 

60 

1 

4 

0.04 

0.017 

0.936 

2.13 

0.1627 

70 

1 

4 

0.03 

0.014 

0.944 

2.13 

0.1622 

80 

1 

4 

0.03 

0.012 

0.951 

2.14 

0.1618 

90 

1 

4 

0.02 

0.011 

0.957 

2.14 

0.1615 

100 

1 

4 

0.02 

0.010 

0.961 

2.14 

0.1613 

110 

1 

4 

0.02 

0.009 

0.964 

2.14 

0.1611 

120 

1 

4 

0.02 

0.008 

0.967 

2.14 

0.1609 

130 

1 

4 

0.02 

0.008 

0.970 

2.15 

0.1608 

66 


TABLE  II  — CONTINUED 


J 

^0 

^1 

^0 

d-P0) 

d-Pi) 

P(E) 

1 

2 

8 

4.32 

0.053 

0.100 

0.44 

0.0761 

2 

2 

8 

2.16 

0.080 

0.238 

0.75 

0.1054 

3 

2 

8 

1.44 

0.144 

0.255 

1.27 

0.1091 

4 

2 

8 

1.08 

0.090 

0.406 

1.20 

0.1147 

5 

2 

8 

0.87 

0.330 

0.202 

2.88 

0.1319 

10 

2 

8 

0.43 

0.181 

0.449 

3.51 

0.0957 

15 

2 

8 

0.29 

0.125 

0.587 

3.76 

0.0910 

20 

2 

8 

0.22 

0.095 

0.670 

3.89 

0.0903 

30 

2 

8 

0.14 

0.064 

0.766 

4.03 

0.0886 

40 

2 

8 

0.11 

0.049 

0.818 

4.10 

0.0849 

50 

2 

8 

0.09 

0.039 

0.852 

4.15 

0.0828 

60 

2 

8 

0.07 

0.033 

0.875  • 

4.18 

0.0815 

70 

2 

8 

0.06 

0.028 

0.892 

4.20 

0.0807 

80 

2 

8 

0.05 

0.025 

0.905 

4.21 

0.0800 

90 

2 

8 

0.05 

0.022 

0.915 

4.23 

0.0796 

100 

2 

8 

0.04 

0.020 

0.923 

4.24 

0.0792 

110 

2 

8 

0.04 

0.018 

0.930 

4.24 

0.0789 

120 

2 

8 

0.04 

0.017 

0.936 

4.25 

0.0786 

130 

2 

8 

0.03 

0.015 

0.940 

4.26 

0.0784 

h 

0.48 

0.83 

1.13 

1.93 

1.88 

4.76 

5.26 
5.54 
5.83 
5.99 
6.09 
6.15 
6.20 
6.24 

6.26 

6.29 

6.30 

6.32 

6.33 


0.0749 

0.0807 

0.0841 

0.0748 

0.0688 

0.0580 

0.0519 

0.0495 

0.0494 

0.0469 

0.0453 

0.0443 

0.0435 

0.0430 


0.0426 


0.0423 


0.0420 


0.0418 


TABLE  II  — CONTINUED 


J 

^0 

U-P0> 

(1-PX) 

io 

P(E) 

1 

3 

15 

7.46 

0.012 

0.018 

0.48 

0.0150 

2 

3 

15 

3.73 

0.066 

0.059 

1.02 

0.0596 

3 

3 

15 

2.49 

0.080 

0.125 

1.37 

0.0305 

4 

3 

15 

1.86 

0.173 

0.112 

2.20 

0.0412 

5 

3 

15 

1.49 

0.122 

0.199 

2.20 

0.0361 

10 

3 

15 

0.75 

0.259 

0.223 

5.22 

0.0371 

15 

3 

15 

0.50 

0.181 

0.368 

5.86 

0.0291 

20 

3 

15 

0.37 

0.139 

0.472 

6.21 

0.0245 

30 

3 

15 

0.25 

0.095 

0.607 

6.60 

0.0203 

40 

3 

15 

0.19 

0.072 

0.687 

6.80 

0.0195 

50 

3 

15 

0.15 

0.058 

0.741 

6.92 

0.0194 

60 

3 

15 

0.12 

0.049 

0.779  ’ 

7.01 

0.0193 

70 

3 

15 

0.11 

0.042 

0.807 

7.07 

0.0184 

80 

3 

15 

0.09 

0.037 

0.829 

7.12 

0.0178 

90 

3 

15 

0.08 

0.033 

0.846 

7.15 

0.0174 

100 

3 

15 

0.07 

0.030 

0.861 

7.18 

0.0171 

110 

3 

15 

0.07 

0.027 

0.873 

7.21 

0.0169 

120 

3 

15 

0.06 

0.025 

0.882 

7.23 

0.0167 

130 

3 

15 

0.06 

0.023 

0.891 

7.25 

0.0165 

TABLE  II— CONTINUED 


a 

J 

^0 

^1 

U-P0) 

(1-PX) 

1 

4 

8 

5.77 

0.215 

0.191 

0.52 

i 

2 

4 

8 

2.89 

0.323 

0.238 

1.10 

r 

3 

4 

8 

1.92 

0.385 

0.255 

1.71 

m 

4 

4 

8 

1.44 

0.264 

0.406 

1.69 

j 

5 

4 

8 

1.15 

0.191 

0.525 

1.61 

10 

4 

8 

0.58 

0.330 

0.449 

4.38 

15 

4 

8 

0.38 

0.234 

0.587 

4.79 

j 

20 

4 

8 

0.29 

0.181 

0.670 

5.01 

30 

4 

8 

0.19 

0.125 

0.766 

5.25 

40 

4 

8 

0.14 

0.095 

0.819 

5.37 

I 

c? 

50 

4 

8 

0.12 

0.077 

0.852 

5.45 

'V 

60 

4 

8 

0.10 

0.064 

0.875 

5.50 

. 

70 

4 

8 

0.08 

0.056 

0.892 

5.54 

| 

80 

4 

8 

0.07 

0.049 

0.904 

5.57 

90 

4 

8 

0.06 

0.043 

0.915 

5.59 

•  • 

100 

4 

8 

0.06 

0.039 

0.923 

5.61 

V  * 

*1 

110 

4 

8 

0.05 

0.036 

0.930 

5.62 

120 

4 

8 

0.05 

0.033 

0.936 

5.63 

130 

4 

8 

0.04 

0.030 

0.940 

5.64 

b 


■j 

72 

5 


P  ( E ) 

0.2031 

0.2620 

0.2461 

0.2361 

0.2325 

0.2329 

0.2249 

0.2242 

0.2136 

0.2097 

0.2078 

0.2067 

0.2060 

0.2055 

0.2052 

0.2049 

0.2047 

0.2045 

0.2044 


TABLE  II — CONTINUED 


!l 

(1-Pq) 

(1-p^ 

^0 

T(E) 

L2 

7.28 

0.051 

0.090 

0.45 

0.0703 

L2 

3.64 

0.143 

0.151 

0.99 

0.1441 

L2 

2.43 

0.151 

0.238 

1.32 

0.1022 

L2 

1.82 

0.264 

0.199 

2.16 

0.1193 

L2 

1.46 

0.191 

0.308 

2.14 

0.1129 

L2 

0.73 

0.330 

0.301 

5.16 

0.1126 

L2 

0.49 

0.234 

0.449 

5.76 

0.0959 

L2 

0.36 

0.181 

0.549 

6.10 

0.0916 

L2 

0.24 

0.125 

0.670 

6.46 

0.0812 

L2 

0.18 

0.095 

0.741 

6.66 

0.0786 

L2 

0.15 

0.077 

0.787 

6.78 

0.0777 

L2 

0.12 

0.064 

0.819 

6.86 

0.0774 

L2 

0.10 

0.056 

0.842 

6.92 

0.0772 

L2 

0.09 

0.049 

0.861 

6.96 

0.0772 

L2 

0.08 

0.043 

0.875 

6.99 

0.0772 

L2 

0.07 

0.039 

0.887 

7.02 

0.0765 

L2 

0.07 

0.036 

0.897 

7.05 

0.0759 

L2 

0.06 

0.033 

0.905 

7.07 

0.0753 

L2 

0.06 

0.030 

0.912 

7.08 

0.0749 

73 


TABLE  II— CONTINUED 


J 

Ho 

Hi 

Ho 

(1-Pq) 

d-Pi) 

Ho. 

P(E) 

1 

5 

25 

12.43 

0.002 

0.003 

0.48 

0.0026 

2 

5 

25 

6.21 

0.014 

0.035 

0.89 

0.0147 

3 

5 

25 

4.14 

0.028 

0.082 

1.24 

0.0107 

4 

5 

25 

3.11 

0.038 

0.130 

1.56 

0.0082 

5 

5 

25 

2.49 

0.080 

0.125 

2.28 

0.0103 

10 

5 

25 

1.24 

0.090 

0.287 

3.58 

0.0085 

15 

5 

25 

0.83 

0.283 

0.189 

8.39 

0.0118 

20 

5 

25 

0.62 

0.221 

0.287 

9.21 

0.0088 

30 

5 

25 

0.41 

0.154 

0.435 

10.15 

0.0062 

40 

5 

25 

0.31 

0.118 

0.535 

10.67 

0.0048 

50 

5 

25 

0.25 

0.095 

0.607 

10.99 

0.0048 

60 

5 

25 

0.21 

0.080 

0.660 

11.22 

0.0040 

70 

5 

25 

0.18 

0.069 

0.700 

11.38 

0.0037 

80 

5 

25 

0.16 

0.061 

0.732 

11.51 

0.0035 

90 

5 

25 

0.14 

0.054 

0.757 

11.60 

0.0034 

100 

5 

25 

0.12 

0.049 

0.779 

11.68 

0.0034 

110 

5 

25 

0.11 

0.044 

0.797 

11.75 

0.0033 

120 

5 

25 

0.10 

0.041 

0.812 

11.80 

0.0033 

130 

5 

25 

0.10 

0.038 

0.825 

11.85 

0.0033 

75 


TABLE  III 


Symmetric  Cutoff  Rate  for  Ideal  Detector 


J 

^0 

V_1 

u-p0> 

(l-p^ 

^0 

P(E)  * 

1 

1 

5 

0.080 

0.124 

0.329 

0.1025 

2 

1 

5 

0.090 

• 

0.287 

0.336 

0.1274 

3 

1 

5 

0.283 

0.189 

0.314 

0.1445 

4 

1 

5 

0.221 

0.287 

0.347 

0.1440 

5 

1 

5 

0.181 

0.368 

0.368 

0.1447 

10 

1 

5 

0.095 

0.607 

0.408 

0.1201 

15 

1 

5 

0.064 

0.717 

0.422 

0.1128 

20 

1 

5 

0.049 

0.779 

0.428 

0.1097 

1 

2 

8 

0.053 

0.100 

0.391 

0.0761 

2 

2 

8 

0.080 

0.238 

0.404 

0.1054 

3 

2 

8 

0.144 

0.255 

0.413 

0.1091 

4 

2 

8 

0.090 

0.406 

0.419 

0.1147 

5 

2 

8 

0.330 

0.202 

0.386 

0.1319 

10 

2 

8 

0.181 

0.449 

0.468 

0.0957 

15 

2 

8 

0.125 

0.587 

0.495 

0.0910 

20 

2 

8 

0.095 

0.670 

0.508 

0.0903 

1 

3 

15 

0.012 

0.018 

0.688 

0.0150 

2 

3 

15 

0.066 

0.059 

0.697 

0.0596 

3 

3 

15 

0.080 

0.125 

0.713 

0.0305 

4 

3 

15 

0.173 

0.112 

0.695 

0.0412 

5 

3 

15 

0.121 

0.199 

0.727 

0.0361 

10 

3 

15 

0.259 

0.223 

0.726 

0.0371 

15 

3 

15 

0.181 

0.368 

0.778 

0.0291 

20 

3 

15 

0.139 

0.472 

0.802 

0.0245 

*  P(E)  from  TABLE  II 


TABLE  III— CONTINUED 


J 

^1 

( 

i-pfl) 

(l-Pj.) 

!o 

P  ( E )  * 

1 

4 

8 

0 

.215 

0.191 

0.149 

0.2031 

2 

4 

8 

0 

.323 

0.238 

0.147 

0.2620 

3 

4 

8 

0 

.385 

0.255 

0.145 

0.2461 

4 

4 

8 

0 

.264 

0.406 

0.160 

0.2361 

5 

4 

8 

0 

.191 

0.525 

0.163 

0.2325 

10 

4 

8 

0 

.330 

0.449 

0.172 

0.2329 

15 

4 

8 

0 

.234 

0.587 

0.189 

0.2249 

20 

4 

8 

0 

.181 

0.670 

0.198 

0.2242 

1 

5 

20 

0 

.014 

0.011 

0.714 

0.0123 

2 

5 

20 

0 

.042 

0.067 

0.732 

0.0434 

3 

5 

20 

0 

.088 

0.101 

0.737 

0.0252 

4 

5 

20 

0 

.132 

0.125 

0.738 

0.0433 

5 

5 

20 

0 

.080 

0.238 

0.752 

0.0339 

10 

5 

20 

0 

.090 

0.406 

0.769 

0.0340 

15 

5 

20 

0 

.283 

0.264 

0.763 

0.0300 

20 

5 

20 

0 

.221 

0.368 

0.800 

0.0243 

*  P(E)  from  TABLE  II 


ft 


.  > . 

TABLE 

IV 

[1 

Avalanche 

Detector , 
P{M0}  = 

Quantized 

P{M1> 

Output 

If 

j 

^1 

N0a 

(1-P0 

(l-P^ 

^0 

P(E) 

■ 

i 

5 

20 

82.36 

0.193 

0.091 

0.59 

0.1419 

2 

5 

20 

30.96 

0.324 

0.121 

1.26 

0.1663 

3 

5 

20 

16.16 

0.376 

0.141 

1.93 

0.1859 

i 

4 

5 

20 

9.77 

0.411 

0.131 

2.67 

0.1388 

5 

5 

20 

6.39 

0.414 

0.130 

3.35 

0.1139 

6 

5 

20 

4.30 

0.418 

0.120 

4.08 

0.1029 

7 

5 

20 

3.13 

0.407 

0.123 

4.70 

0.0739 

a 

8 

5 

20 

3.01 

0.367 

0.161 

4.99 

0.0771 

9 

5 

20 

5.62 

0.283 

0.286 

4.48 

0.0808 

S 

V 

».  • 

10 

5 

20 

7.51 

0.229 

0.393 

4.09 

0.1052 

>  t 

i 

c* 

1 

5 

25 

101.87 

0.141 

0.059 

0.59 

0.0997 

»/! 

2 

5 

25 

39.11 

0.265 

0.094 

1.25 

0.1252 

v" 

3 

5 

25 

20.84 

0.333 

0.106 

1.95 

0.1447 

‘‘,1 

| 

4 

5 

25 

12.77 

0.361 

0.110 

2.64 

0.1003 

5 

5 

25 

8.50 

0.370 

0.112 

3.32 

0.0829 

v* 

6 

5 

25 

5.93 

0.387 

0.094 

4.13 

0.0692 

7 

5 

25 

4.18 

0.372 

0.103 

4.71 

0.0497 

8 

5 

25 

3.17 

0.367 

0.101 

5.38 

0.0357 

9 

5 

25 

2.65 

0.374 

0.090 

6.17 

0.0289 

•  .' 

10 

5 

25 

4.87 

0.280 

0.205 

5.46 

0.0352 

11 

5 

25 

6.68 

0.224 

0.315 

4.91 

0.0530 

* 

12 

5 

25 

8.09 

0.187 

0.408 

4.49 

0.0605 

I 


TABLE  IV— CONTINUED 


N0a 

U-P0> 

U-Pi) 

P(E) 

150.36 

0.125 

0.066 

0.56 

0.0957 

61.39 

0.231 

0.104 

1.19 

0.1254 

34.12 

0.306 

0.125 

1.86 

0.1331 

21.83 

0.359 

0.130 

2.57 

0.1100 

15.09 

0.379 

0.143 

3.21 

0.1120 

10.98 

0.410 

0.129 

4.01 

0.1104 

8.27 

0.406 

0.138 

4.62 

0.0812 

6.39 

0.414 

0.130 

5.36 

0.0668 

4.98 

0.438 

0.104 

6.32 

0.0505 

3.90 

0.443 

0.095 

7.12 

0.0435 

3.21 

0.413 

0.119 

7.46 

0.0342 

2.79 

0.431 

o.ioo  • 

8.43 

0.0258 

3.39 

0.363 

0.166 

8.03 

0.0330 

5.05 

0.289 

0.276 

7.10 

0.0441 

6.39 

0.256 

0.337 

6.81 

0.0484 

7.51 

0.229 

0.393 

6.54 

0.0521 

8.48 

0.207 

0.442 

6.30 

0.0588 

80 


TABLE  IV- -CONTINUED 


J 

^0 

*1 

z 

o 

fit 

(1“Po) 

(1-P2) 

1 0 

P  ( E ) 

1 

8 

40 

185.04 

0.075 

0.040 

0.55 

0.0575 

2 

8 

40 

76.23 

0.177 

0.072 

1.19 

0.0851 

3 

8 

40 

43.20 

0.253 

0.092 

1.86 

0.0918 

4 

8 

40 

27.81 

0.307 

0.099 

2.57 

0.0705 

5 

8 

40 

19.44 

0.334 

0.110 

3.24 

0.0712 

6 

8 

40 

14.31 

0.350 

0.116 

3.90 

0.0707 

7 

8 

40 

10.88 

0.370 

0.108 

4.67 

0.0507 

8 

8 

40 

8.50 

0.370 

0.112 

5.31 

0.0424 

9 

8 

40 

6.76 

0.380 

0.102 

6.10 

0.0372 

10 

8 

40 

5.43 

0.376 

0.104 

6.74 

0.0268 

11 

8 

40 

4.40 

0.377 

0.099 

7.47 

0.0189 

12 

8 

40 

3.82 

0.386 

0.086 

8.34 

0.0136 

13 

8 

40 

3.10 

0.363 

0.105 

8.67 

0.0120 

14 

8 

40 

2.75 

0.383 

0.084 

9.73 

0.0084 

15 

8 

40 

3.42 

0.324 

0.142 

9.23 

0.0099 

16 

8 

40 

4.87 

0.280 

0.205 

8.74 

0.0122 

17 

8 

40 

6.06 

0.231 

0.303 

7.78 

0.0196 

18 

8 

40 

7.07 

0.208 

0.356 

7.45 

0.0223 

19 

8 

40 

7.93 

0.198 

0.375 

7.57 

0.0215 

81 


TABLE  IV— CONTINUED 


J 

vo 

^1 

^Oa 

d-P0) 

(1-PX) 

P(E) 

1 

10 

40 

195.01 

0.102 

0.052 

0.56 

0.0769 

2 

10 

40 

82.36 

0.193 

0.091 

1.17 

0.1053 

3 

10 

40 

47.26 

0.266 

0.113 

1.83 

0.1047 

4 

10 

40 

30.98 

0.324 

0.121 

2.53 

0.0887 

5 

10 

40 

21.83 

0.359 

0.130 

3.22 

0.0941 

6 

10 

40 

16.16 

0.376 

0.141 

3.86 

0.0939 

7 

10 

40 

12.41 

0.394 

0.139 

4. 57 

0.0759 

8 

10 

40 

9.77 

0.411 

0.131 

5.34 

0.0662 

9 

10 

40 

7.84 

0.420 

0.125 

6.09 

0.0624 

10 

10 

40 

6.39 

0.414 

0.130 

6.70 

0.0487 

11 

10 

40 

5.23 

0.413 

0.129 

7.37 

0.0394 

12 

10 

40 

4.30 

0.418 

0.120 

8.15 

0.0338 

13 

10 

40 

3.80 

0.430 

0.105 

9.08 

0.0268 

14 

10 

40 

3.13 

0.407 

0.123 

9.40 

0.0209 

15 

10 

40 

2.79 

0.431 

0.100 

10.54 

0.0151 

16 

10 

40 

3.01 

0.367 

0.161 

9.98 

0.0203 

17 

10 

40 

4.43 

0.320 

0.225 

9.44 

0.0212 

18 

10 

40 

5.62 

0.283 

0.286 

8.97 

0.0284 

19 

10 

40 

6.63 

0.253 

0.342 

8.55 

0.0302 

20 

10 

40 

7.51 

0.229 

0.393 

8.18 

0.0366 

TABLE  V 


Symmetric  Cutoff  Rate  for  Avalanche  Detector 

P{Mq}  =  P{M1> 


J 

^0 

<1-P0> 

d-P^ 

fo 

P  ( E )  * 

1 

5 

20 

0.193 

0.091 

0.243 

0.1419 

2 

5 

20 

b.324 

0.121 

0.259 

0.1663 

3 

5 

20 

0.376 

0.141 

0.280 

0.1859 

4 

5 

20 

0.411 

0.131 

0.333 

0.1388 

5 

5 

20 

0.414 

0.130 

0.399 

0.1139 

6 

5 

20 

0.418 

0.120 

0.478 

0.1029 

7 

5 

20 

0.407 

0.123 

0.546 

0.0739 

8 

5 

20 

0.367 

0.161 

0.578 

0.0771 

9 

5 

20 

0.283 

0.286 

0.518 

0.0808 

10 

5 

20 

0.229 

0.393 

0.464 

0.1052 

1 

5 

25 

0.141 

0.059 

0.332 

0.0997 

2 

5 

25 

0.265 

0.094 

0.352 

0.1252 

3 

5 

25 

0.333 

0.106 

0.382 

0.1447 

4 

5 

25 

0.361 

0.110 

0.435 

0.1003 

5 

5 

25 

0.370 

0.112 

0.499 

0.0829 

6 

5 

25 

0.387 

0.094 

0.586 

0.0692 

7 

5 

25 

0.372 

0.103 

0.648 

0.0497 

8 

5 

25 

0.367 

0.101 

0.714 

0.0357 

9 

5 

25 

0.374 

0.090 

0.775 

0.0289 

10 

5 

25 

0.280 

0.205 

0.724 

0.0352 

11 

5 

25 

0.224 

0.315 

0.661 

0.0530 

12 

5 

25 

0.187 

0.408 

0.603 

0.0605 

*  P(E)  from  TABLE  IV 
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TABLE  V 


Symmetric  Cutoff  Rate  for  Avalanche  Detector 


—  * 

P<M0> 

=  P(MX} 

J 

^0 

<1-P0> 

(l-P^ 

** 

fo. 

P (E )  * 

1 

5 

20 

0.193 

0.091 

0.243 

0.1419 

2 

5 

20 

t) .  324 

0.121 

0.259 

0.1663 

3 

5 

20 

0.376 

0.141 

0.280 

0.1859 

4 

5 

20 

0.411 

0.131 

0.333 

0.1388 

5 

5 

20 

0.414 

0.130 

0.399 

0.1139 

6 

5 

20 

0.418 

0.120 

0.478 

0.1029 

7 

5 

20 

0.407 

0.123 

0.546 

0.0739 

8 

5 

20 

0.367 

0.161 

0.578 

0.0771 

9 

5 

20 

0.283 

0.286 

0.518 

0.0808 

10 

5 

20 

0.229 

0.393 

0.464 

0.1052 

1 

5 

25 

0.141 

0.059 

0.332 

0.0997 

2 

5 

25 

0.265 

0.094 

0.352 

0.1252 

3 

5 

25 

0.333 

0.106 

0.382 

0.1447 

4 

5 

25 

0.361 

0.110 

0.435 

0.1003 

5 

5 

25 

0.370 

0.112 

0.499 

0.0829 

6 

5 

25 

0.387 

0.094 

0.586 

0.0692 

7 

5 

25 

0.372 

0.103 

0.648 

0.0497 

8 

5 

25 

0.367 

0.101 

0.714 

0.0357 

9 

5 

25 

0.374 

0.090 

0.775 

0.0289 

10 

5 

25 

0.280 

0.205 

0.724 

0.0352 

11 

5 

25 

0.224 

0.315 

0.661 

0.0530 

12 

5 

25 

0.187 

0.408 

0.603 

0.0605 
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J 

^0 

(1-P0) 

(1-Pi) 

^0_ 

P(E)  * 

1 

6 

12 

0.375 

0.229 

0.062 

0.3019 

2 

6 

12 

0.476 

0.240 

0.063 

0.3244 

3 

6 

12 

0.528 

0.228 

0.072 

0.3369 

4 

6 

12 

0.541 

0.220 

0.093 

0.2943 

5 

6 

12 

0.583 

0.173 

0.130 

0.2588 

6 

6 

12 

0.572 

0.180 

0.156 

0.2434 

7 

6 

12 

0.396 

0.377 

0.127 

0.2643 

8 

6 

12 

0.316 

0.488 

0.112 

0.2812 

1 

6 

18 

0.247 

0.131 

0.170 

0.1895 

2 

6 

18 

0.362 

0.166 

0.180 

0.2181 

3 

6 

18 

0.420 

0.182 

0.188 

0.2338 

4 

6 

18 

0.462 

0.165 

0.223 

0.1937 

5 

6 

18 

0.469 

0.161 

0.271 

0.1685 

6 

6 

18 

0.476 

0.149 

0.331 

0.1549 

7 

6 

18 

0.465 

0.153 

0.385 

0.1254 

8 

6 

18 

0.470 

0.144 

0.442 

0.1020 

9 

6 

18 

0.327 

0.323 

0.363 

0.1323 

10 

6 

18 

0.266 

0.428 

0.320 

0.1577 

*  P(E)  from  TABLE  IV 


TABLE  V— CONTINUED 


J 

^0 

^1 

u-p0> 

(l-P^ 

^0 

P(E)  * 

1 

8 

32 

0.125 

0.066 

0.338 

0.0957 

2 

8 

32 

0.231 

0.104 

0.373 

0.1254 

3 

8 

32 

0.306 

0.125 

0.381 

0.1331 

4 

8 

32 

0.359 

0.130 

0.400 

0.1100 

5 

8 

32 

0.379 

0.143 

0.422 

0.1120 

6 

8 

32 

0.410 

0.129 

0.470 

0.1104 

7 

8 

32 

0.406 

0.138 

0.513 

0.0812 

8 

8 

32 

0.414 

0.130 

0.570 

0.0668 

9 

8 

32 

0.438 

0.104 

0.643 

0.0505 

10 

8 

32 

0.443 

0.095 

0.702 

0.0435 

11 

8 

32 

0.413 

0.119 

0.726 

0.0342 

12 

8 

32 

0.431 

0.100 

0.777 

0.0258 

13 

8 

32 

0.363 

0.166 

0.764 

0.0330 

14 

8 

32 

0.289 

0.276 

0.701 

0.0441 

15 

8 

32 

0.256 

0.337 

0.673 

0.0484 

16 

8 

32 

0.229 

0.393 

0.645 

0.0521 

17 

8 

32 

0.207 

0.442 

0.620 

0.0588 

*  P(E)  from  TABLE  IV 
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TABLE 
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J 

^1 

<l-p0)  (1-Pl) 

*0 

P(£)  * 

1 

8 

40 

0.075 

0.040 

0.453 

0.0575 

2 

8 

40 

0.177 

0.072 

0.493 

0.0851 

3 

8 

40 

0.253 

0.092 

0.507 

0.0918 

4 

8 

40 

0.307 

0.099 

0.530 

0.0705 

5 

8 

40 

0.334 

0.110 

0.555 

0.0712 

6 

8 

40 

0.350 

0.116 

0.590 

0.0707 

7 

8 

40 

0.370 

0.108 

0.640 

0.0507 

8 

8 

40 

0.370 

0.112 

0.684 

0.0424 

9 

8 

40 

0.380 

0.102 

0.739 

0.0372 

10 

8 

40 

0.376 

0.104 

0.780 

0.0268 

11 

8 

40 

0.377 

0.099 

0.823 

0.0189 

12 

8 

40 

0.386 

0.086 

0.866 

0.0136 

13 

8 

40 

0.363 

0.105 

0.880 

0.0120 

14 

8 

40 

0.383 

0.084 

0.913 

0.0084 

15 

8 

40 

0.324 

0.142 

0.903 

0.0099 

16 

8 

40 

0.280 

0.205 

0.885 

0.0122 

17 

8 

40 

0.231 

0.303 

0.833 

0.0196 

18 

8 

40 

0.208 

0.356 

0.809 

0.0223 

19 

8 

40 

0.198 

0.375 

0.817 

0.0215 

*  P(E) 

from 

TABLE 

IV 
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J 

^o 

^1 

U-P0) 

a-px) 

fo. 

P(E)  * 

1 

10 

30 

0.179 

0.104 

0.241 

0.1414 

2 

10 

30 

0.273 

0.143 

0.278 

0.1703 

3 

10 

30 

0.345 

0.160 

0.286 

0.1721 

4 

10 

30 

0.388 

0.177 

0.287 

0.1564 

5 

10 

30 

0.420 

0.182 

0.297 

0.1644 

6 

10 

30 

0.452 

0.168 

0.328 

0.1614 

7 

10 

30 

0.471 

0.158 

0.365 

0.1329 

8 

10 

30 

0.482 

0.149 

0.410 

0.1128 

9 

10 

30 

0.475 

0.153 

0.452 

0.1081 

10 

10 

30 

0.476 

0.149 

0.499 

0.0920 

11 

10 

30 

0.485 

0.136 

0.553 

0.0731 

12 

10 

30 

0.456 

0.161 

0.576 

0.0759 

13 

10 

30 

0.478 

0.137 

0.630 

0.0619 

14 

10 

30 

0.372 

0.257 

0.569 

0.0741 

15 

10 

30 

0.327 

0.323 

0.540 

0.0772 

16 

10 

30 

0.291 

0.382 

0.514 

0.0873 

17 

10 

30 

0.262 

0.435 

0.486 

0.0975 

18 

10 

30 

0.237 

0.482 

0.462 

0.1009 

♦  P(E) 

from 

TABLE 
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TABLE  V— CONTINUED 


J 

is 

il 

(1-P0) 

d-Pi) 

is. 

P  ( E )  * 

1 

10 

40 

0.102 

0.052 

0.389 

0.0769 

2 

10 

40 

0.193 

0.091 

0.438 

0.1053 

3 

10 

40 

0.266 

0.113 

0.451 

0.1047 

4 

10 

40 

0.324 

0.121 

0.463 

0.0887 

5 

10 

40 

0.359 

0.130 

0.477 

0.0941 

6 

10 

40 

0.376 

0.141 

0.495 

0.0939 

7 

10 

40 

0.394 

0.139 

0.530 

0.0759 

8 

10 

40 

0.411 

0.131 

0.573 

0.0662 

9 

10 

40 

0.420 

0.125 

0.619 

0.0624 

10 

10 

40 

0.414 

0.130 

0.659 

0.0487 

11 

10 

40 

0.413 

0.129 

0.701 

0.0394 

12 

10 

40 

0.418 

0.120 

0.749 

0.0338 

13 

10 

40 

0.430 

0.105 

0.796 

0.0268 

14 

10 

40 

0.407 

0.123 

0.815 

0.0209 

15 

10 

40 

0.431 

0.100 

0.854 

0.0151 

16 

10 

40 

0.367 

0.161 

0.841 

0.0203 

17 

10 

40 

0.320 

0.225 

0.816 

0.0212 

18 

10 

40 

0.283 

0.286 

0.789 

0.0284 

19 

10 

40 

0.253 

0.342 

0.762 

0.0302 

20 

10 

40 

0.229 

0.393 

0.734 

0.0366 

*  P(E)  from  TABLE  IV 
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